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: LIMITING VALUES OF GRAVITATIONAL AND 
MAGNETIC ANOMALIES DUE TO A SUB- 
TERRANEAN STRUCTURE BOUNDED 
BY A SINGLE DIFFERENTIAL 
SURFACE* 


J. W. FISHERT 


ABSTRACT 


Calculations are made of gravitational and magnetic anomalies supposed due to 
local variations in the form of a single differential surface separating an overlying rock 
from an underlying one, the effective thickness of the latter being infinite. Both layers 
are homogeneous in density or magnetic susceptibility and only those magnetic effects 
are considered which are due to induction in the earth’s magnetic field. 

The maximum possible values of these anomalies and of their first and second hori- 
zontal derivatives are calculated under these simplified conditions, and it is pointed 
out that the numerical readings of these quantities, provided by survey results, may 
lead to useful estimates of the extreme depth and proportions of the structure respon- 


sible. 





A fundamental problem of geophysics is that of deducing from 
surface measurements of the gravitational and magnetic fields at 
various points over a given region and kind of underground rock, 
structures which may account for the observed variations from nor- 
mal. One of the simplest assumptions we might be permitted to make 
in certain cases is that the observed surface anomalies may be account- 
ed for by projections of an underlying basement rock into an overly- 
ing one, each of the two layers being assumed homogeneous in density 
and magnetic susceptibility and the lower of the two supposed to 
extend downwards to an infinite depth. Even with such simplifying 
assumptions, the problem of determining the form of the boundary 
surface separating the two layers, so as to fit in with the observed sur- 
face fields, may be very difficult and tedious to solve, depending as it 


* Manuscript received July 1940. 
+ Geophysical Department, Anglo-Iranian Oil Co., Ltd. 
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does to a large extent on trial and error. In a previous paper,' an ex- 
perimental method of computation has been described which might 
assist in this process. However, in certain practical cases, all that may 
be required of the survey is an estimate of the depth of the upper sur- 
face of the basement layer and the vertical extent of the folds in its 
bounding surface. Now it can be shown that, whatever the nature of 
the folding of the differential surface separating the two layers .{ 
rock, there will be certain maximum possible values for the anomalies 
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and their horizontal derivatives on the surface of the ground, and that 
these limiting values will depend solely on the assumed difference in 
density or susceptibility, and on the greatest and least depths of the 
bounding surface below the surface of the ground. A knowledge of 
these theoretical maxima enables us alternatively, by comparison 
with the maximum observed values of the corresponding quantities, 
to fix certain limits to the depth of the bounding surface. An upper 
limit, for instance, for the shallowest depth can be found from the 
maximum observed value of the second and higher horizontal deriva- 
tives of the vertical gravity anomaly or from the first and higher hori- 
zontal derivatives of the vertical magnetic anomaly. The geometrical 
form of those structures of the bounding surface which would lead to 
these theoretical maxima are often widely different than any which 
might be expected as approximations to actual geological structures, 
but the existence of such absolute limits is none the less important on 
that account. These limiting values will now be calculated in order, 
starting with the gravitational anomaly. 


1 J. W. Fisher: Gropuysics, Vol. 5, No. 1, Jan. 1940. 
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The density excess of the basement layer over that of the surface 
layer will be denoted by p, assumed to be positive. The corresponding 
relative magnetic susceptibility will be denoted by k. The structure of 
the boundary surface is supposed confined between two horizontal 
planes, the upper of which is at a depth / and the lower at a depth h’ 
below the surface, which is taken as the xy plane, the z axis being sup- 
posed vertically downwards (Fig. 1.). 


A. VERTICAL GRAVITATIONAL ANOMALY 


Consider a layer of vertical thickness Az at a depth z. An element 
of matter of area dxdy at a point (x,y,z) of this layer will, if occupied 
by basement rock, contribute to the vertical gravitational force at a 
surface point P (£,n,o) an amount 


GpzAzdxdy/[z? + (« — £2? + (y — n)?}8, 


where G is the gravitational constant. Since z is always positive, this 
contribution must always be positive. The contribution of the whole 
layer at depth z to the vertical gravitational anomaly at the surface 
point P is 


As ) C A ff zdxdy ( ) 
Pi z T 5 : . 
$ ¥ [32 + (vw — £2 + (y — 0)? )8”? 


and the maximum effect at P will be obtained by filling up all the 
space between z=/ and z=h’ by the basement rock. This gives the 
maximum possible value of the anomaly dg, viz., 


(6g2)max = 27pG(h’ — h). (2) 





If (6g-)obs is the maximum anomaly actually observed, we may be 
sure that 


ian es 


(2a) 
2mpG. ‘ 


B. FIRST HORIZONTAL DERIVATIVE OF THE 
VERTICAL GRAVITATIONAL ANOMALY 


By differentiating (1) under the integral sign we have 


, s(x — £)dxdy 
— A(ég.) = 30 
de (dg2) = 3 ff [s? + (2 — 8 + (9 — (3) 
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where o=GpAz. If we put «/=x—£, y’=y—y7, the contribution of the 
layer at depth z becomes 


7) A(s =f f ° x'dx'dy’ 
at gz) [22 oe x’? .’ yy’? |5/2 


To find, say, the greatest positive value which the layer can contrib- 
ute, we include all elements for which the integrand is positive. Thus 
if x’ is positive, we can fill up the whole of the positive half plane for 
any depth z. This leads to a step-like structure with its edge at x’=o, 
extending to + © in the ydirection and of height h’—A. Integrating 
over the above range of variables, we find for the maximum first 
horizontal derivative of the vertical gravitational anomaly the value 


S = 20G log h’/h. (4) 





This value applies to points along the y’ axis immediately above 
the edge of the step. In this case, therefore, the optimum structure 
reduces to a two-dimensional one. If |.Sobs| is the greatest observed 
absolute value of the first horizontal derivative, we arrive at the limit 

h’ a obs | 
log — 2 
h 2pG 





(4a) 





C. SECOND HORIZONTAL DERIVATIVE OF THE 
VERTICAL GRAVITATIONAL ANOMALY 


Differentiating (3) with respect to & we get 


0? ge if ig ~~ e 
hae = ixdy 
ae? (6g) = 30 Ke =a ree yan 


fas (y x + 2?) dx'dy’. (5) 
ae 4 vy? oh. 3? |7/2 


The sign of the integrand is that of 4x%”—(y’-+-2*), and the regions of 
positive and masiaan contributions are separated by the hyperbola 
4x"—y?—2=0, or y?=4x"—2, with asymptotes y’= + 2x’ (Fig. 2). 
For a given z the elements in the shaded zone will provide a negative 
contribution. We integrate (5) first with respect to x’ between limits 
x’=o0 and x/=}(2+y”)!? and then with respect to y’ from o to o, 
This gives the negative contribution of one shaded quadrant, and 


therefore } of the total effect of the layer thus bounded. We find 




















VALUES OF GRAVITATIONAL AND MAGNETIC ANOMALIES 5 


12y1/2 yy"? op 52 _ 4x’? 


1 @ x b(z*+y 
— — — A(ég,) = sos f ay’ f dx 
4 o0& 0 0 [et fe of? an 





Finally, integrating over z between the limits / to h’, we derive for the 
maximum absolute value of the second horizontal derivative, 


48rpG / 1 I I I 
E = ——= (= = =) = 2.7pG (— —— ~). (6) 
25V5\h h’ h h’ 


Both positive and negative values of any derivative have of course the 
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same absolute numerical value, since the effect of a completely filled 
up layer is zero at any surface point. If now | T obs| is the greatest ab- 
solute value observed for the second horizontal derivative, we can say 
that 





ieee | ? I | T obs | 
—, ie, —>- 
h h’ 2.7pG h 2.7pG 
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This result is useful because it leads to an upper limit to /, the shallow- 
est depth of the structure, viz., 


h < 2.70G/ | T ss . (7) 


It is of interest to compare this result with that which would apply to 
an infinitely long ridge with straight sides, since such a structure would 
be a more likely approximation to those geological formations which 
may be met with in practice. Treating this problem as a two-dimen- 
sional one, but adopting the same method as described above, we find 
that the maximum effect is produced by a truncated ridge, the sides of 
which are planes with their line of intersection on the surface and 
inclined 60° to the horizontal. The numerical factor in this case turns 
out to be 2.598 as compared with 2.7 in the case of the three-dimen- 
sional structure calculated above. 





D. HORIZONTAL DIRECTING TENDENCY (DIFFERENTIAL 
CURVATURE) RECORDED BY THE 
EOTVOS TORSION BALANCE 


Apart from the horizontal gradient, the limiting value of which 
has been calculated under B, the Torsion Balance provides values of 
the H.D.T., which is the product of ‘g’ and the difference between the 
principal curvatures of the equipotential surface passing through the 
point of measurement. If &, » denote axes directed along those of 
principal curvature at the observation point P, we define the H.D.T. 


by 
= =) 
dn ae J’ 


where U denotes the gravitational potential. The quantity 
aU 
Onde 





vanishes in this case. 
Our problem, therefore, is to calculate the greatest absolute value 


of 
e *U) 
On? 0 ‘ 


Introducing as before «’=x—£, y’=y—7 we find that a layer at depth 
z contributes an amount 


U 


x/2 — y* : 
30 - —— dx'dy’. 
[«’2 + y+ g? |5/2 
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The integrand is positive over the shaded area bounded by 
y’=+2’ (Fig. 3). Integrating over this whole region gives 


) x’ 4/2 == 5 i 
r20 f av’ f 19 1" 215/2 dy’ 
0 0 [x 2 oe y 2 + g? |8/2 


, rt z'*dg! 4pGAz 
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A final integration over z between the limits /, h’ gives for the maxi- 
mum absolute value of the H.D.T. the value 


H.D.T.max = 40G log h'/h. (8) 


FIG. 3 


It is interesting to note that a similar calculation in the case of a two- 
dimensional structure shows a maximum absolute value of half the 
above amount i.e., 2pG log h’/h. This two-dimensional structure is in 
the form of a long truncated ridge with sides sloping at 45°, and inter- 
secting in the surface, the top and bottom surfaces being at depths 
hand h’ respectively below the surface. 

It would of course be possible to calculate limits for the third and 
higher horizontal derivatives of 5g. and so obtain upper limits for 
h?, h® etc., from the observed maxima of these higher derivatives. In 
practice, however, it would not be feasible to obtain from gravity 
traverses sufficiently accurate values of derivatives higher than the 
second. Using the same method we proceed to calculate the limiting 
values of the vertical magnetic anomaly and its horizontal deriva- 
tives. 








8 J. W. FISHER 


E. VERTICAL MAGNETIC ANOMALY 


Let k denote the relative susceptibility of the basement magnetic 
rock, supposed magnetized inductively by the earth’s field H. Let H, 
denote the vertical component of H and H, the horizontal com- 
ponent, supposed directed along the x axis. If Aw represents an element 
of solid angle subtended at a surface point P(é, 7, o) by an element of 
rock in the horizontal layer of vertical thickness Az at a depth 2z, the 
contribution of this element to the vertical anomaly at P will be 


0 0 
k (1.j.20 + H.2.d0) 


Writing 

Aw = zdxdy / [22 + (x — &)? + (y — 0)?” 
and integrating over the whole horizontal layer at depth z, we get for 
its contribution :— 


222 — x’? — vy’? + 302%! 
kH Az dq dx'dy’ 
(32? + x2 + y’?)5/2 ’ 


with a=H,/H, and y’=y—n, x’=x—£& as before. The boundary 
between positive and negative contributions will be the circle 
x?+y"— 3azx’— 22?=0, and the greatest positive effect of the layer 
will therefore be found by integrating over the interior of this circle, 
at all points of which the integrand is positive. Had we considered the 
vertical component only of the earth’s field, the section of optimum 
structure would have been a circle of.radius z\/2 with its center ver- 
tically beneath the point of maximum vertical anomaly. The optimum 
structure would then be a truncated cone with apex at the surface 
point. This case, which admits of immediate integration over x’, y’ 
and z, leads to a maximum anomaly of 


4nvkH, 
sV3° 


The effect of the horizontal component of the earth’s field is to dis- 
place the center of the circular cross-section in the positive x direction 








log h'/h = 2.42kH, log h'/h. 


: que ; , , 3a\? 
through a distance —— and to increase its radius to z 2+( 3 
2 


2 
The integration over x’, y’ is best performed graphically after trans- 
forming to polar co-ordinates x’=r cos ¢, y’=r sin ¢. This allows im- 
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mediate integration over 7 for a fixed ¢. The lower limits of the r 
integrals in this case are zero and the upper limits are found graphi- 
cally for different values of ¢. The area of the curve obtained by plot- 
ting @ and the definite integral gives the total required effect of the 
horizontal layer. If we take a= H,/H.=0.42, the layer at depth z con- 
tributes 2.57 kH./z and so the maximum vertical anomaly for a range 
of z between / and h’ would be 


(6H z)max = 2.57kH,z log h'/h. (9) 


As in the case of the first horizontal derivative of the gravitational 
anomaly, this can be used to provide an upper limit to the ratio h’/h. 


F. FIRST HORIZONTAL DERIVATIVE OF THE VERTICAL MAGNETIC 
ANOMALY PERPENDICULAR TO MAGNETIC MERIDIAN 
(MERIDIAN PARALLEL TO x OR TO & AXIS) 


We find 


d (=) fia ff Sy at = 5" = FP ee? os 
. ae a x’dy’. 
On\ kHz h (s2 + a/2 4 y/2)7/2 . 


Thus, for a positive value of the derivative, we require the sign of 
42°—x'?— y+ 5azx’ to be the same as that of y’. The numerator of 
the integrand is, as in the previous cases, homogeneous in powers of 
x’, y’, z, and the final integral over z will be of the form 


W’ dz I I 
sf Sd 
h 2° h h’ 


where A is the numerical value of the double integral in the x’y’ 
plane when z has the value 1. The center of the boundary circle 
x+y’ — sax’ —4=0 in the plane at unit depth (z=1) is displaced by 
a distance 1.05 in the x’ direction, if a is given the value 0.42 (Fig. 4). 
The region of integration is the interior of the semicircle for positive 
y’, and for negative y’ it will be over the external region to infinity. 
The contribution to A of the interior of the semicircle is found graphi- 
cally to be 2.35. The effect of the complete half plane for y’ positive is 
found by immediate integration to be 2. Thus the positive y’ region 
outside the semicircle must contribute —o.35. The contribution of 
the corresponding region for y’ negative will therefore be +0.35, and 
the total value of A for the shaded area is accordingly 2.35+0.35= 
2.70. The maximum value of the first horizontal derivative perpendic- 
ular to the magnetic meridian will therefore be 2.7 RH(1/h—1/h’). 
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As in the case of the second horizontal derivative of the gravitational 
anomaly, this gives an upper limit to hd. 


PARALLEL TO MAGNETIC MERIDIAN 


For 0/0¢(5H./kH.) positive, the cross-section for unit depth (z= 1) 
has the form shown in Fig. 5, the shaded regions being occupied by 
magnetic matter. The equation to the boundary between positive and 


y 


/ 


\ 











FIG. 4 


negative regions is now 32” [4—(x?+y”) |—3a[1+y"—4x”|=0. On 
carrying out the integrations graphically, the extreme value is found 
to be the same as in the former case viz., 2.70 kH.(1/h—1/h’). 

To summarize, if the magnetic meridian is along the x axis and 
the ratio of horizontal to vertical component of the earth’s field is 
0.42, then 


|= ox) | - |< 0x | ee ee ) 
0g ‘ ™ On ' max ——o {+ = ne 


This value may be compared with the value which would be found for 
the vertical part of the earth’s field acting alone. By immediate inte- 
gration this is found to be 2.58 kH.(1/h—1/h’). 








Gi 
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G. SECOND HORIZONTAL DERIVATIVE OF THE 
VERTICAL MAGNETIC ANOMALY 


The limiting value of 0?/d£(6H.) either parallel or perpendicular to 
the magnetic meridian could be found by the same general method, 
but in order to avoid very laborious calculations, the value is calcu- 
lated in this case taking account of the vertical component of the 











earth’s field only. Inclusion of the horizontal component should in- 
crease this value somewhat, just as it does in the case of the anomaly 
itself and its first derivative. 

Under these conditions (a=o), we have 


02 





6H. 
oe OH) 


a’ g4— 2722 y/2+ 222 y’2+4 3x +4x/4— y’4 
= 31H, [ as { {+ ; - : A : a dx'dy’. 
‘ [« 24 y 2 +4. 3? |9/2 





For regions where the numerator of the integrand is positive, we have 
contributions to the negative value of the second derivative. This 
leads to a complicated type of cross-section of the form shown in 
Fig. 6, taking z unity. 

The shaded areas indicate those portions of the horizontal cross- 
section at unit depth which contribute to the negative second deriva- 
tive at a surface point immediately over the origin. The double integral 
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over the central shaded area has the numerical value 5.81, and that 
over the outer shaded area has the value 0.176. Thus the effect of the 
central part of the structure is about 30 times that of the outer part. 
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A layer of thickness Az at a depth z will therefore contribute an 
amount 5.99 kH.(Az/z*) to the negative second derivative. Integrating 
over z we Have 


| : or) | ene : (11) 
oO aa 7 


or, if T denotes the maximum observed absolute value of the second 
derivative, h<1/3kH./T, and this may provide a useful indication 
of an upper limit for the shallowest depth of the structure. The same 
numerical factor applies to the positive values of the second deriva- 
tive, so that only the absolute values need be taken into account. 
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THE GULF GRAVIMETER* 
R. D. WYCKOFF 


ABSTRACT 


The design features of the Gulf gravimeter! are presented accompanied by a brief 
discussion of general gravimeter design problems. In addition to the description of the 
unastatized instrument which is the type finally selected for commercial development, 
brief consideration is given to two astatized instruments including a bifilar gravimeter 
the construction of which was carried to completion. 


INTRODUCTION 


Gravimetric surveying has, until recent years, suffered from a 
lack of adequate instrumentation properly to adapt it to the role for 
which it is best fitted. To alleviate this situation arising from the 
fact that the torsion balance is inherently a detailing device, ill 
adapted to reconnaissance, consideration was given to direct gravity 
measurements whereby the limitations of the gradient method could 
be eliminated. This led to the early development of the pendulum 
method and while commercial use of the apparatus has been rather 
limited in the industry at large, a considerable amount of such work 
was done within the Gulf Oil Corporation organization. However, 
while this work demonstrated that creditable performance may be 
achieved with properly designed pendulum equipment, it was recog- 
nized at the start that greater speed and the attendant decrease in 
cost would be necessary if gravity surveying was to maintain its de- 
served scope of operation. Acknowledging this fact and recognizing 
that inherent limitations in the pendulum method precluded attain- 
ment of this end, early in November 1932 Dr. Eckhardt instigated in- 
vestigation of the gravimeter design problem at the Gulf laboratory. 
Early in 1933 definite plans had been outlined and laboratory work 
was well underway. 

With the view to minimize delay in providing a field-worthy in- 
strument the plan involved the simultaneous investigation and con- 
struction of two distinctly different models well known at the time 
and now commonly referred to as the astatized and unastatized types. 
By following this procedure the necessity of making a final choice as 


* Presented at the Annual Meeting, Chicago, April 9, 1940. 
+ Gulf Research & Development Co., Pittsburgh, Pa. 
1 A. Hoyt, U.S. Patent No. 2, 131,737, issued 10/4/38, appl. 8/5/35. 
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to type was postponed until a rational selection could be based upon 
actual experience with each type. At the end of 1934, after two 
years’ elapsed time, an unequivocal selection was possible and in 
March 1935 a successful field trial of one of the instruments was made 
in the coastal area of Texas. This instrument, supplemented by im- 
provements where experience indicated the need, represents the 
gravimeter which has been used by the Gulf organization since the 
first eight were put into operation and completely supplanted the 
pendulum equipment in 1936. 

-It is the purpose of this discussion to present the results of the 
above investigation. Attention will be confined primarily to a specific 
instrument but in so doing it will be necessary to consider briefly the 
instrumental requirements and the salient features of the two types 
mentioned above. No attempt will be made to provide a comprehen- 
sive discussion of gravimeter design problems and an exhaustive com- 
parison of the various instruments now in use would be quite unjusti- 
fied. 


GENERAL CONSIDERATIONS IN GRAVIMETER DESIGN 


SENSITIVITY REQUIREMENTS 


Experience amply demonstrates that even for reconnaissance 
prospecting purposes, a satisfactory gravimetric survey requires an 
instrument having a probable error not exceeding o.1 milligals and 
preferably less. The lowest limit to be assigned is for obvious reasons 
determined by practical considerations involving the accuracy with 
which elevation and topographic corrections? can be made and by 
erratic surface and near-surface influences arising from density varia- 
tions. While the magnitude of these spurious and uncontrollable effects 
will vary depending upon the geological and terrain environment, it 
is probably safe to say that little practical significance can be attached 
to the attainment of an accuracy exceeding 0.05 mg. unless such ac- 
curacy is available without sacrifice of operational convenience. Ex- 
cept for very special applications the utilitarian advantage of an 
instrument of such precision lies in the assurance that the normal ob- 
servational errors will be below the “noise level’’ due to terrain and 
topographic uncertainties, and that from normal probability consid- 
erations, the number of stations exceeding. 0.1 mg. error will be 
correspondingly reduced. 


2 Elevation corrections:—free air correction=0.094 mg./ft. or including Bouguer 
correction=0.07 mg./ft. approx. 
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Accepting for the purpose of discussion that a probable error not 
exceeding o.1 mg. is desired, and recalling that in the gravimeter-type 
instrument the total force of gravity, 980,000 mg., is operating, the 
minimum sensitivity required is approximately 10~" g. Since this re- 
quirement of measuring a force to better than one part in 10,000,000 
must be accomplished by a portable instrument easily manipulated 
in the field under any adverse terrain and weather conditions en- 
countered, and must be sufficiently rugged to withstand reasonably 
rough handling and transportation, some of the difficulties in design 
will be evident. In fact absolute measurements of any force with that 
accuracy have not and are not likely soon to, be attained. 

The gravimeter is retrieved from this apparently hopeless situation 
by the use of proper operating methods since unpredictable instru- 
mental drift is readily determined by the simple expedient of return- 
ing at suitable time intervals to reobserve a previously occupied sta- 
tion. If the “drift curves’ obtained in this manner are relatively 
smooth and contain no sudden breaks, reliable observations may be 
possible even though random long-period fluctuations or drift may be 
several orders of magnitude greater than the probable error specified. 
Evidently then the term gravimeter is somewhat of a misnomer. A 
strictly applicable term would be gravity comparator for this designa- 
tion eliminates the implication of long-term stability usually asso- 
ciated with indicating instruments and its use might eradicate the 
skepticism of the uninitiated when a geophysicist speaks glibly of the 
accuracy of gravity measurements in terms of parts in one hundred 
million. Here we will at least qualify the term gravimeter with the 
understanding that the instrument is a comparator. 


UNASTATIZED TYPE (LINEAR SYSTEM) 


In our original investigations starting in 1933 all efforts were soon 
concentrated upon instruments of the type involving a constant mass 
the weight of which is balanced by a comparison force derived from 
a suitable spring. In its simplest form such a balance comprises merely 
a mass supported by a helical spring in conjunction with a device to 
indicate the equilibrium position of the mass. 

It is evident that if in an elemental device of the above type (Fig. 
1a), S is the total displacement of the spring-mass system under the 
action of total gravity g, then the sensitivity of the device as a gravim- 
eter will be: 


AS = S-Ag/g. (x) 
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That is, a change in the gravitational force of 1077 g is accompanied 
by a change in the position of the weight amounting to 10~” times the 
total initial displacement. Hence, if we assume 20 cm. (5 inches) as 
the maximum displacement of a spring which can be incorporated 
conveniently in a portable instrument we obtain a displacement: 


AS = 2X 10-§cm. when Ag = 1077 g. 


This displacement amounts to about 1/30 of the fringe separation 
in an interferometer or about 100 times smaller than the limiting reso- 
lution of a precision micrometer microscope. Since we have been speak- 
ing of minimum sensitivity requirements it is evident that without 
some amplification of the motion of the suspended mass, neither the 
microscope nor interferometer is applicable. In fact, for other reasons 
we may eliminate the latter device from any further consideration. 

Adopting the micrometer microscope as a suitable measuring 
device because of its convenience, we may assume that a displace- 
ment of 0.001” or 2.5X10-* cm. is directly observable with ease. But 
we have seen that in a portable gravimeter of the elemental spring- 
balance type a displacement of 2X10 cm. must be detected with 
certainty, so it is evident that an amplification of at least 1000 
must be incorporated between the mass system and the observed 
index point. While a portion of this amplification may be obtained by 
optical means, full advantage of optical levers can be secured only 
when a rotational motion is available. This is a disconcerting situation 
which is met in some instruments’ by mechanical conversion systems 
which amplify directly by levers or convert the linear motion of the 
suspended mass into a rotational motion of small mirrors. Since any 
auxiliary attachments of this kind introduce hinge or pivot points 
and ligaments contributing extraneous restoring forces, and neces- 
sarily introduce multiple constraint in the moving system, tending 
to aggravate level sensitivity, a strenuous effort was made to retain 
the simplicity of the simple spring balance. The problem was finally 
solved in the Gulf instrument by incorporating the rotational feature 
within the spring itself as will be explained later. 


ASTATIZED TYPE (NON-LINEAR SYSTEM) 


Thus far we have mentioned only the method of applying amplifi- 
cation between the moving mass and the reading device and noted that 
for adequate sensitivity an amplification X1000 or greater must be 
obtained, 
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An alternative method involves increasing the actual motion of 
the moving mass by the method commonly called astatization. If it 
were possible to introduce an absolutely constant restoring force to 
balance say 99% of the total weight of the moving system then the 
apparent sensitivity of the balance operating as a gravimeter could 
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Fic. 1. Schematic diagram of unastatized (a) and astatized (b) spring balances. 
K=spring constant. The schematic force/displacement curves (b) illustrate the prin- 
ciple of astatization. 


be increased 100 fold. While this ideal astatization of a spring-balance 
system has never been attained, it has been simulated by innumerable 
well known methods exemplified by the Galitzin vertical seismograph 
and its numerous modifications or the bifilar gravimeter.* 

Fig. 1 shows schematically a comparison of the unastatized spring 
balance and a typical astatized system though it should be emphasized 
that these figures are intended merely to indicate the essential differ- 
ence in the mechanics of the two types. 


3 That is, Tomaschek and Schaffernicht, ‘“The Temporal Variation in the Force of 
Gravity,” Annalen d. Physik, Vol. 15, p. 787, (1932). 
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In the case of the elemental balance of Fig. 1a, the forces acting 
are obviously the constant weight Mg supported by the opposing 
spring tension which increases linearly with displacement at a rate 
proportional to the spring constant K. As indicated by the Figure the 
equilibrium displacement is shown by the intersection of the curves 
KS and Mg. It is also evident from the Figure that the sensitivity of 
the device as a gravimeter will increase as angle 6 decreases; a result 
which, because of the linearity of curve KS, can be accomplished only 
by increasing the equilibrium elongation of the spring. Obviously this 
is merely a restatement of Eq. 1 but is intended to emphasize the fact 
that in any spring balance system having a strictly linear restoring- 
force curve, KS, the initial or equilibrium displacement uniquely 
determines its gravimetric sensitivity.* 

On the other hand, while the ideally astatized system would be 
linear, practically attainable systems are of the non-linear type indi- 
cated in Fig. rb, the turning moment or torque functions being de- 
pendent upon the geometry of the system. Thus, by attachment of 
the spring to the lever L’ rather than on the center line Z between the 
mass and knife-edge, the instrument represented by Fig. rb may be 
adjusted so that its equilibrium point is anywhere along the restoring- 
torque curve: 


T = K[L’ cos (6 + 4) |[L’ sin (6 + ¢) + So]. (approx.) (2) 


Thus angle 8 may be made arbitrarily small and the sensitivity arbi- 
trarily great, up to the point of instability when B=o. Hence by de- 
vices equivalent to Fig. rb it is possible to simulate the behavior of 
extremely long springs and thereby attain very high sensitivities. 
However, it is evident that with each change in adjustment, the sensi- 
tivity and hence the scale value of such an instrument will be changed. 
Moreover, even for a fixed adjustment, the non-linear feature pre- 
cludes a strictly linear scale, although over a limited range which may 
be extended by the use of further geometrical complexity and auxiliary 
spring arrangements, reasonable linearity may be achieved. This in- 
herent non-linearity of the so-called astatized devices is inconvenient 
and while certain design features may alleviate this handicap, it con- 


4 Evidently the actual physical elongation of a spring may be made somewhat less 
than the effective elongation by the expedient of pre-stressing the spring when forming 
it. However, this merely suppresses its zero and in speaking of its displacement we refer 
to effective elongation. 
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stitutes the principal objection to this type. Moreover, it practically 
dictates the design of a null rather than a deflection type instrument. 

It is of interest to note that the gravimetric sensitivity of any 
spring balance is determinable from the period of the system, and this 
provides a direct and simple criterion whereby a comparison of various 
types may be made. Thus in any simple harmonic system such as a 
weight in equilibrium with a linear restoring force the period, 7, in 
seconds is given by 


T = 20(M/Kg)"2, (3) 


where: 
M = Mass in gms. 
K =Spring constant in gms/cm. 
g= Acceleration of gravity in dynes. 
But M/K=Sb, the initial or equilibrium displacement of the spring, 


therefore 
T = 2n(So/g)!/?, (4) 


T? = (40°/g)So. (5) 


Hence from Equations 1 and 5, we may say that the gravimetric 
sensitivity of any spring balance is proportional to the square of its 
period and moreover this applies not only to the linear but to the asta- 
tized systems’ which cannot be compared on the basis of any other 
simple criterion. 

Consider for example a simple spring balance, Fig. 1a, having a 
total spring displacement of 20 cm. and a sensitivity, mentioned 
earlier, of 2X10-® cm. per 10~“g. The period of the system is (Eq. 4) 
approximately o.g seconds. Now in order to obtain an increase in 
sensitivity of 100 fold, the device of Fig. rb may be adjusted to a 
period of g seconds. With such an adjustment the instrument simu- 
lates the performance of an elemental balance having an equilibrium 
displacement of more than 2000 cm. and yet in the astatized system 
the actual elongation of the spring need not exceed a few centimeters. 


5 This is true, of course, only to the extent that we may neglect inertial components 
in the system other than those upon which gravitational force is acting. For example, 
the rotational moments of the beam system of Fig. 1b represent extraneous inertial 
factors which increase the period of the system without affecting the gravimetric sensi- 
tivity. Obviously the period may be corrected for these spurious effects before computing 
the sensitivity. 
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Fic. 2. Galitzin Type Gravimeter. 
. Knife-edge lift and clamp control 
. Main-spring tension adjuster 
. Main-spring 
. Mass 
. Fixed mirror for multi-reflection 
6. Moving mirror (top of knife-edge) 
7. Knife edge 
8. Lifter and clamp arm 
g. Entrance optical slit and exit window to photo recorder 
10. Optical system for level (15) 
11. Micro-adjustment of spring angle 
12. Main spring support 
13. Main spring support hinge 
14. Auxiliary-spring micrometer control 
15. Pendulum type level 
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This method of effectively amplifying the displacement of the moving 
-mass, together with the use of a very moderate optical amplification, 
results in a compact device of high sensitivity and presents a most 
alluring picture the disadvantages of which, along with the advan- 
tages, can only be fully appreciated through experience. 


GRAVIMETER CONSTRUCTION 


ASTATIZED TYPE 


Mention was made that our plan of investigation involved con- 
struction of both astatized and unastatized types. In view of existing 
information regarding the former type and the immediate need for 
an instrument of high sensitivity in which to test the performance of 
various springs and spring materials, an experimental model was con- 
structed which was nothing more than a modified Galitzin vertical 
seismograph with a knife-edge support in place of the usual hinge. 
It was adapted to continuous photographic recording to facilitate 
examination of the performance characteristics of the springs. The 
design is that shown schematically in Fig. rb with necessary refine- 
ments to permit micro-adjustments of spring tension and alignment 
angle of the spring with respect to the axis of the moving system. The 
experimental instrument is shown in the photograph, Fig. 2. 

With a properly thermostated container which maintained the 
temperature constant to 0.001°C. this instrument easily could be 
adjusted to a period of 10 or 12 seconds and periods in excess of 25 
seconds were maintained for 24 hours or more. In Fig. 3 the circled 
points represent deflections as read from a photographic record after 
subtracting a linear drift and indicate the excellent performance ob- 
tained when the period was adjusted to approximately 22 seconds. 
These data are presented merely to show the extremely high sensi- 
tivity obtainable through the use of the astatization principle. 

While this instrument was provided with an arresting and clamping 
device, for various and rather obvious reasons no serious attempt was 
made to adapt it to field use. It served admirably, however, as a 
spring tester and incidentally to develop tangible experience and famil- 
iarity with the problems involved in instruments intended to operate 
at sensitivities in excess of 10~"g. 

Fig. 4 shows the important features of a bifilar gravimeter designed 
as a portable instrument and having adequate sensitivity for gravi- 
metric surveys. It represents our attempt to build a field-worthy in- 
strument incorporating the astatization principles. While this bifilar 
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Fic. 3. Gravity variations due to tidal forces and associated phenomena presented 
to show the high sensitivity obtainable with compact instruments of the astatized 


type. 


design and construction was carried through to completion it was 
never given a field test although laboratory operation indicates reason- 
ably satisfactory performance and ruggedness. Of special interest is 
the calibrating weight and the auxiliary sylphon-bellows arrangement 
to provide pressure variations for the buoyancy method of calibration. 
These self-contained calibrating devices were provided in order con- 
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. 4. Bifilar Gravimeter. 
. Clamp control 
. Spring-tension supports most of mass, torsion opposes bifilar torque 


Bifilar suspension—Drawn in to make visible 


. Clamp arm 
. Mirror on moving system 


Torsion head adjuster 
Calibrating weight 


8. Mass 


9. 
10. 


Ligament operating clamp-arm 
Damping magnet 


veniently to determine the scale value of the instrument for any 
desired adjustment and to detect the drift in sensitivity which might 
easily occur in a device of this type. It appeared to be the logical 


method 


of combating the inconvenient non-linear characteristics of 


this type of instrument. 
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UNASTATIZED TYPE 


We have already discussed the sensitivity requirements of an ele- 
mental spring balance and concluded that a magnification of approxi- 
mately X1000 would be required between the moving mass and the 





Fic. 5. Complete moving system of Gulf Gravimeter. 


observed index point. The desirability of eliminating any lever systems 
or other mechanical amplifying schemes will also be apparent. 

Fig. 5 shows the spring-mass system of the Gulf gravimeter and 
it will be noted that it retains the extreme simplicity and desirable 
features of the elemental spring balance of Fig. 1a—a helical spring 
supporting a moving mass with no ligaments, levers, or other attach- 
ments. The unique feature is, of course, the use of a ribbon helix 
which, in the spring illustrated, produces a rotation of about 1 second 
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of arc per 10-7 change in g. This is equivalent to an initial deflection 
under total g of approximately 2880° or 8 revolutions. The rotation 
of the suspended mass is observed by means of an attached mirror 
which in conjunction with a special multiple-reflection mirror arrange- 
ment permits reading the straight-through, or undeflected beam 
together with multiple reflections readable beyond the 6th. Using the 
4th reflection and an optical lever 22” in length, o.1 milligal corre- 
sponds to a 0.001” displacement of the index under the micrometer- 
microscope. Experience indicates that the instrument can be read 
to about ten times this accuracy although, of course, the usable sen- 
sitivity is somewhat less. 

Certain very important operating advantages are available with 
the simple moving system of Fig. 5, for in addition to adequate sensi- 
tivity when used in conjunction with moderate optical amplification, 
the adjustment of the instrument has no marked effect on the scale 
value or calibration constant. The total range of the micrometer read- 
ing device is 25 to 30 mg. depending upon the particular spring in- 
volved but it is a very simple matter to adjust the instrument for 
operation in any area regardless of the magnitude of g, and the cali- 
bration constant is unaffected by such adjustments. The moving sys- 
tem is suspended from a torsion head provided with micrometer ad- 
justment which is available to the operator at any time by merely 
inserting a special adjusting shaft through an opening in the top con- 
trol panel. Thus, if in the course of a day’s work gravity variations 
exceed the range of the scale the operator merely shifts the zero of the 
instrument by means of the torsion head adjustment and proceeds 
without appreciable loss of time or change in calibration constant. 
Instruments which have been in service for several years have checked 
their original calibration within 0.5% and better. 

Another convenient feature of the moving system arising however 
from the nature of the spring material rather than its design, is the 
fact that the rate of drift is adjustable within reasonable limits. While 
the principal criterion of a good spring is a uniform rate of creep, it is 
nevertheless very convenient to work with an instrument having a 
low rate of drift. In the Gulf gravimeter the direction of creep is a 
function of the stress in the spring and by merely adjusting the mass to 
the proper value so as to operate at or near the point of reversal, the 
creep rate may be made relatively small. A typical instrument oper- 
ating in a fixed location and observed by the author for some 30 days, 
showed a drift of about 0.3 mg./24 hrs., although, of course, in actual 
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field operation the drift would depend somewhat upon operating 
conditions. 
DAMPING 

To facilitate accurate reading and rapid operation the moving sys- 
tem must be properly damped. This is necessary not only to eliminate 
natural oscillation of the system immediately after it is unclamped 
but to prevent the building up of sustained oscillations due to micro- 
seisms, wind, traffic, and other disturbances. In the Gulf instrument 
the circumstances are unusually favorable for their elimination. Or- 
dinary short-period seismic disturbances cause vertical displacements 
of the weight and are not observable in the optical system since the 
index point is deflected by rotation and only inappreciably by trans- 
lational movements of the mirror. Now the longitudinal period of the 
system in Fig. 5 is about 0.7 to 1 second whereas the rotational period 
is approximately 7 seconds. This disparity in natural period between 
the vertical and rotational modes, of which only the latter affects the 
optical system, results in a very inefficient conversion of short period 
seismic disturbance to scale deflection. Of course long-period earth- 
quake and similar disturbances are observed but due to the small 
coupling between vertical and rotational modes at high frequencies 
the ordinary short-period disturbances unless of great magnitude are 
effectively filtered out by the moving system itself. 

Electromagnetic damping is used and is obtained in the conven- 
tional manner by a metallic vane moving through the air gap of a 
permanent magnet. The damping vanes indicated in Fig. 5 are ex- 
tended to a large radius with the view to increase both the damping 
and the rotational moment of inertia of the system. It is evident that 
damping is applied to both vertical and rotational motions of the sys- 
tem but most effectively to the important rotational mode. 


LEVEL SENSITIVITY 


One of the most disconcerting problems encountered in all gravim- 
eters is the elimination of level sensitivity. This is particularly true 
of certain types which inherently are as sensitive to level as to changes 
in gravity. In the Gulf instrument it is evident that the absence of 
any mechanism attached to the moving system, that is the absence 
of any multiple restraint on the moving mass, eliminates the most 
potent sources of level sensitivity. However, in its action as a plumb- 
bob there is bending of the spring immediately adjacent to the support 
and this bending is partially converted into rotation by the ribbon 
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helix. By the simple expedient of tapering the helix at this point the 
rotational tendencies are effectively compensated. Alternatively, other 
supporting means which permit bending unaccompanied by rotation 
have been used. 

We are accustomed to speak of level sensitivity in terms of level 
tolerance, that is the amount the instrument may be shifted off level 
in the direction of maximum sensitivity before the instrument deflec- 
tion reaches a value of o.1 mg. This level tolerance differs in individual 
instruments ranging from 10 to 12 seconds of arc in the worst case to 
more than 60 seconds, and is due partly to uncompensated effects 


<< ‘| 


Fic. 6. Schematic diagram of multiple-reflection optical system in Gulf Gravimeter. 
Multiple images of slit are observed by micrometer- microscope. 


SLIT MICROMETER 
b MICROSCOPE 





PARTIAL REFL. MIRRORS 


incident to the above mentioned bending of the spring near the point 
of support but primarily to spurious effects in the damping system. 

The instrument is provided with a pair of 10” spirit levels with the 
usual go° azimuthal arrangement. Since this permits leveling the in- 
strument to an accuracy of a few seconds of arc no difficulty is en- 
countered in normal operations in view of the instrumental level 
tolerance stated above. The level adjustment is accomplished by 
means of the conventional three-screw arrangement the only novel 
feature being in the construction of these screws whereby the threads 
as well as a ball thrust bearing operate in an oil bath and dust is per- 
manently eliminated. 

OPTICAL SYSTEM 


The arrangement of the optical system is shown schematically in 
Fig. 6. The moving and fixed mirrors are combined with the lens sys- 
tem so as to eliminate unnecessary loss of light. The flat sides of the 
plano-convex lenses are good optical flats and are coated by an evapo- 
ration process with a partially reflecting aluminum film which, for 
optimum light intensity at the eyepiece, has a transmission coefficient 
of approximately 12%. 

An important feature of the optical system is the fact that shifting 
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Fic. 7. Gulf Gravimeter. 


1. Thermometer 11. Electromagnetic damping 
2. Levels 12. Torsion-head adjuster 
3. Temperature compensator 13. Optical slit 
4. Level-sensitivity minimizer 14. Alignment adjustment for mass 
5. Ribbon spring 15. Mass 
6. Clamp-arm stop 16. Moving mirror 
7. Clamp 17. Fixed mirror 
8. Ligament operating clamp-arm 18. Prisms 
9. Clamp-arm 19. Clamp-control gear reduction 
10. Clamp shaft (3) 


of the light source (slit) or the eyepiece cannot introduce errors since 
the zero position, or undeflected image, is always available as a refer- 
ence or true-zero index. A micrometer reading from this index to any 
one of the multiple-reflection images provides a measure of the rota- 
tional deflection of the moving system which is free from any zero- 
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Fic. 8. Gulf Gravimeter control panel. 


1. Clamp control 4. Micrometer-microscope 


. Light source for slit 5. Level windows 
3. Opening for torsion-head control 6. Opening for pressure equalizing control 


N 


shift errors. Moreover, since reflections beyond the 6th order are 
readable and the displacement of the image is two times the order of 
the reflection, a reading on the 4th may be independently checked by 
a simultaneous reading on say the 6th reflection and thus elimination 
of gross reading errors is possible. The micrometer-microscope used 
is a conventional type reading to 0.oo01”/division. The average Gulf 
instrument has a scale constant of 100 divisions/mg. on the 4th reflec- 
tion so that no difficulty is experienced in reading to better than 0.05 
mg. under routine conditions. 
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CLAMPING DEVICE 


Fig. 7 shows the complete assembly of the instrument proper in 
which the arrangement of the clamping mechanism may be seen. The 
three V-grooved clamping lugs on the moving mass (Fig. 5) are en- 
gaged by three ball-tipped pins on the clamping arms (Fig. 7). The 
clamp-arms are actuated by a single control-shaft through reduction 
gearing which operates the three clamp-shafts and thence through 





Fic. 9. Top view of Gulf Gravimeter assembly—Control panel and top insulation 
disc removed. 


1. Thermostat well 5. Exit window—indicator light 
2. Heater shell 6. Entrance slit—indicator light 
3. Clamp control 7. Torsion head adjuster 


4. Pressure equalizing valve 


ribbon ligaments attached to the ends of the clamp-arms which pivot 
on hinges. Adjustments provided are adequate to permit precision 
adjustment of the clamps and if proper care has been taken the moving 
system may be clamped without seriously disturbing the images in 
the reading microscope. 

Fig. 8 shows the top plate or operating panel of the instrument. It 
will be noted that the clamping controi is provided with a safety de- 
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vice which causes a lug to project beyond the edge of the top plate 
when the moving system is unclamped. In this position it is impos- 
sible to place the dust cover on the instrument and if the operator has 
neglected to clamp the system he is warned before the instrument 
is moved. 

TEMPERATURE CONTROL 


While the material used in the spring has a low temperature coeffi- 
cient and is supplemented by a temperature compensator, no attempt 
has been made to operate without temperature control. Fig. 9 shows 
a top view of the instrument with the upper insulation pad and top 
plate removed. 

The instrument proper is suspended within a thermostatically con- 
trolled heater shell which in turn is surrounded by insulating material 
of low thermal conductivity protected on the outside by a thin alumi- 
num cylinder. The heater shell is arranged to surround the instrument 
completely and the heater winding is distributed over the entire shell 
including the top and bottom of the cylinder. Temperature control 
is obtained by means of a conventional mercury-type thermostat 
mounted in a well attached to the outside of the heater can. A suitable 
relay operating with a few milliamperes control-current from the 
thermostat serves to handle the heater current supplied by a 12-volt 
battery. The instrument is, of course, maintained at a temperature 
higher than the maximum ambient temperature and the total heat 
dissipation of the system is about 0.4 watt/°C. temperature differ- 
ence. This moderate power drain permits the use of a small battery 
when unusual portability is required, the total weight of such auxiliary 
temperature control equipment including battery is 45 lbs. as com- 
pared with 85 lbs. for normal use. 

Tests have shown the control to be good to .oo1°C. for short pe- 
riods or when not transported and 0.02°C. over long periods in field 
operation. 

MAGNETIC SENSITIVITY 


In the early Gulf instruments the moving system exhibited a slight 
permanent magnetic moment resulting in a spurious sensitivity to 
horizontal magnetic fields. In the earth’s field this sensitivity was such 
as to require orientation of the instrument to within +5 or ro degrees 
in a specified azimuth to be within a o.1 mg. tolerance. Since it is a 
relatively simple matter to orient within less than 1 degree with an 
ordinary compass obviously no serious errors were introduced from 
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this cause. However, in more recent instruments this inconvenience 
has been eliminated, the magnetic sensitivity in any ordinary magnetic 
fields being entirely negligible. 

FIELD OPERATION 


In common with most gravimeters, the speed of operation is very 





Fic. 10. Gulf Gravimeter showing yoke near center of gravity for supporting it in 
transportation, dust cover, leveling screws, carrying handles and temperature control 
cable. 

Total weight of instrument proper 25 lbs. 

Total weight of complete instrument 95 lbs. 

Total weight of battery and control equipment 45-85 lbs. 


rapid, the time involved in actually reading the instrument being 
negligible in comparison with the time required to transport it be- 
tween stations. For reasons of policy we have never attempted to 
observe from a tripod -lowered from the truck bed so at every station 
the instrument is removed from the truck and set on three small stakes 
driven into the ground to serve as a convenient support. This proce- 
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dure adds some to the labor involved but is a negligible factor in the 
over-all speed of operation. As illustrative of the speeds which are ob- 
tainable, a single instrument with operator and one helper has ob- 
served over 100 stations with 1 mile spacing in an eight-hour period 
although, of course, normal operation would be less than half that 
number. 

Since 1936 more than 200,000 stations have been observed by the 
Gulf gravimeters and statistics on these routine data indicate that the 
probable error of a single observation is a little less than 0.05 mg. For 
example, a set of data obtained in Kuwait, Arabia, where the operator 
repeated certain stations twenty or thirty times during several months, 
shows a probable error of a single observation of 0.04 mg. and the 
distribution follows a typical probability curve. These operations 
were under rather adverse conditions. Other examples taken at ran- 
dom from routine field data but under more favorable circumstances 
show for instance that repeat observations involving two different 
instruments on a group of 15 or 20 stations resulted in a maximum 
difference of 0.1 mg., a mean deviation of 0.04 mg., and a probable 
error of a single observation of 0.02 mg. A second group of stations 
showed a maximum difference between observations of 0.08 mg., a 
mean deviation of 0.03 mg., and a probable error of 0.02 mg. 

It is of interest to note that the accuracy of gravimeters is sufficient 
to permit their use in obtaining gradient values comparable in ac- 
curacy to the torsion balance. Thus an array of three gravimeter sta- 
tions forming a right angle with legs only 300 feet in length if observed 
with a readily obtained prevision of 0.02 mg. will provide a gradient 
measurement equal in accuracy to that obtainable with the torsion 
balance having a probable error of 4e. For certain special problems we 
have used the gravimeter in this manner with good success. 
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THE GULF UNDERWATER GRAVIMETER* 
T. B. PEPPERt 


ABSTRACT 


A brief description of the Gulf underwater gravimeter is presented showing the 
salient design features whereby a regular Gulf gravimeter is incorporated in a device 
permitting precision gravity measurements in coastal waters several hundred feet in 
depth. 

The instrument is adapted to be lowered to the bottom on a cable and through re- 
mote control devices all necessary adjustments are made from the surface. The reading 
is obtained photographically on a continuously moving film which permits accurate 
observations even when moderate seismic disturbances cause oscillation of the index 
point. 


In the last decade, as existing oil-producing areas have encroached 
steadily upon the water bodies of the Gulf Coast and othez territories, 
the necessity of prospecting in areas covered by water has been keenly 
felt by most of the oil companies. As a result various prospecting meth- 
ods have been tried with varying amounts of success. Seismograph 
operations under such conditions proved to be successful from the start 
but rather laborious and slow for reconnaissance work. 

Pendulums and gravimeters were not particularly adaptable in 
areas covered by water because of the nature of the instruments which 
require a relatively firm support, and particularly pendulums whose 
accuracy is very seriously impaired by excessive seismic disturbance. 
The Gulf gravimeter, while not affected by the adverse conditions as 
seriously as pendulums, was nevertheless affected to the extent that 
it was difficult to read unless mounted on a rather firm foundation or 
tripod. 

By a system of two tripods, one for the operator to stand on and 
one for the gravimeter, it was possible to make observations on calm 
days in water up to eighteen feet in depth. But when the water was 
rough it was necessary to observe for a few hours to get a satisfactory 
reading if, indeed, any was obtained. When obtained, however, the 
accuracy of the observations was better than o.1 milligal as proven 
by repeat observations. The task of setting these tripods, each of which 
weighed 1500 lbs was done by derricks on barges and was slow, making 
the entire operation quite costly. 

It was evident that operations could be simplified and conse- 


* Presented at the Annual Meeting, Chicago, IIl., April 9, 1940. 
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quently expedited if an instrument could be designed to operate on 
the bottom instead of on tripods and that the savings effected by in- 
creased speed of operations would justify a comparatively large ex- 
penditure for such an instrument. 

Many schemes for attaining such an end were suggested and con- 
sidered. Some of the most promising ones are listed: 


(1) Bell or housing containing the operator and instrument. 

(2) Submergeable barge unit containing only the instrument. 

(3) Pressure housing type; instrument to be leveled directly by 
gravity; all other controls to be either automatic or manipu- 
lated by remote control from the surface. 

(4) Same as (3) except that the instrument is to be leveled by re- 
mote control through the use of level indicators. 


Also combinations of any of the above, as well as one which had 
the leveling mechanism outside of the pressure housing. 

In judging the schemes just mentioned, cost, simplicity, safety, 
ease and speed of operation, weight, and durability were kept in 
mind. 

The first scheme, utilizing a diving bell containing operator and 
instrument was rejected because it had none of the qualities desired. 
The only reason it was considered at all was that on first thought it 
appears to be the most direct solution until further consideration 
shows clearly the details and hazards involved. 

The second scheme, utilizing a submergeable barge containing the 
instrument only, eliminated the hazards of using a diving bell and 
obviated the necessity of unloading and loading on board a boat, but 
added the complication of submerging and raising the barge. However, 
it did not present much promise for speed, and had all of the problems 
of schemes (3) or (4) which involve a pressure housing with remote 
controls. 

The third scheme, involving an instrument leveled directly by 
gravity and enclosed in a pressure housing solved the problem of 
leveling, but since it was considered necessary to know with certainty 
whether the instrument was level or not, it gave promise of being as 
complicated as one leveled by remote control. Additional complica- 
tions were involved in damping oscillations of the instrument in the 
pressure housing and clamping it while moving from one location to 
the next. Also it would involve considerable weight because of the 
space necessary for the damping mechanism. If a liquid was used for 
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damping, there would always be the possibility of getting it on the 
optical system and other delicate parts. If electro-magnetic damping 
was used, there would be the possibility of stray magnetism affecting 
the operation of the instrument since parts of it are affected by strong 
fields. 

The fourth scheme, involving an instrument enclosed in a pressure 
housing and leveled by remote control, had possibilities of being the 





Fic. 1. Gulf Underwater Gravimeter sealed with all 
control equipment necessary to operate. 


lightest of all systems excepting the one with the leveling mechanism 
outside of the pressure housing. Protection afforded the leveling mech- 
anism when placed inside of the pressure housing was considered an 
advantage which more than offset the disadvantage of added weight. 
Since this arrangement came closest to satisfying all requirements, 
it was chosen as the type to be built. 

Before an instrument could be designed along these lines a device 
was required which would indicate the level position to within less 
than 10 seconds of arc (see Fig. 2). A very durable combination was 
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devised which utilized a level vial in combination with photo-voltaic 
cells with no amplification. In this combination one cell is used on 
each end of the vial to compensate for light-source voltage fluctua- 


























Fic. 2. This sketch shows the relation between level indicators and the axes about 
which the instrument is revolved to level it. Vial B’ is independent of the state of level 
on A’ and will not be affected if the level of A’ is changed. A’ will be affected if B’ is 
changed to greater or lesser degree dependent on how much the instrument is tilted. 


tions. The cells are connected in parallel but in such manner as appar- 
ently to short each other, i.e., positive to negative, and leads are run 
from the junctions to a 25-0-25 micro amp. meter. Illumination for 
both ends of the vial is from a common source. When the meter indi- 
cates zero, the instrument is level. At the zero position on the micro- 
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ammeter, no appreciable change is noted when the light source voltage 
is varied from zero to any voltage the lamp is able to withstand. This 
system has a sensitivity of about 3 seconds of arc per micro-ampere 
on the most sensitive portion of the indicator response-curve with 
about —10 to +10 micro-amps. total deflection when 8 volts is im- 





Fic. 3. Gulf Underwater Gravimeter with the top half of the pressure housing re- 
moved showing the top-plate assembly and leveling mechanism. Size of available electric 
motors and other essentials almost proved to be a limiting factor on the size of the unit. 


pressed on the lamp. The unit indicates the direction off level even 
when it is out of level more than 40°; a feature of considerable impor- 
tance. 

Leveling is accomplished by a motor-driven gimbal arrangement 
in which the two components of leveling motion are sufficiently inde- 
pendent to permit leveling both directions simultaneously. For this 
purpose a compact arrangement is used in which the maximum lever- 
ages are obtained in a minimum space (see Fig. 3). B and C are light 
sources for the level units. D, E, F, and G are photo-voltaic cells which 
are part of the level indicator system. L is the motor drive for one com- 
ponent of level and M is its associated leveling unit. 
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The maximum available leveling motion is 15° off center in any 
direction. To tilt from 15° in one direction to 15° in any other direction 
requires less than three minutes including time for all minor adjust- 
ments. 

In order to appreciate the precision with which the leveling motors 
have to be stopped it need only be noted that each revolution of the 
motors represents three seconds of arc, and that on high speed they 
run at approximately 15,000 RPM. On fine adjustments they must 
operate at less than 10 RPM. The motors used for this purpose are 
permanent magnet miniature-train motors which while admirably 
suitable in other respects unfortunately have low starting torque, a 
characteristic which complicated matters and necessitated the design 
of a motor driven contactor device with variable duration of contact. 
Electrical impulses from these contactors jerk the motors for a small 
and definite part of a revolution per contact closure. With this sort 
of a device it is possible to level the instrument while making a record 
without destroying the accuracy of observation. 

The level mechanism is so constructed that it has less than one 
second back lash, i.e., when the controls are reversed, there is less than 
3 turn of the motors before the level starts changing. If the instru- 
ment is tilted more than the leveling mechanism can handle the 
motors move the leveling mechanism until 15° off center is reached. 
At this limit a contact is made with a circular contact ring thus oper- 
ating a relay turning off the level motors. In order to make the motors 
operable again it is necessary to press a button (center top on control 
panel) continuously and run the motors in the reverse direction to 
bring it again within the 15° circle. 

Clamping the moving system of the gravimeter is accomplished 
by a permanent-magnet motor, J, working through a gear box, K. 
Warning signals are actuated by a contact on the clamp shaft which 
operates a relay on the operator’s control panel. The relay turns on a 
red pilot lamp (center bottom of the control panel), when the instru- 
ment is unclamped and a green one (center top) only when it is 
clamped completely. 

Recording is accomplished by photographic means on 35 mm. film. 
It is of the continuous type and provides a complete record of the 
instrument deflection from the time the moving system is unclamped 
until it is again clamped. Accuracy of recording in spite of possible 
film shrinkage in processing is obtained by fiducial lines photographed 
simultaneously with the images used for recording gravity. These 
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fiducial lines are put on the record from the back side of the film 
while the gravity indications are photographed from the front. Film 
is drawn through the camera by a permanent-magnet motor, I, driv- 
ing through a 27,000 to 1 worm gear train, H, which builds up suff- 
cient friction to stop the motor when the end of the film is reached. 
An ammeter in this circuit indicates an increase in current at the con- 
trol panel, informing the operator that the film is all used. The camera, 
A, which is a remodelled Leica, can then be removed and reloaded 
when the instrument is brought on board the boat. Film moves 
through the camera at about 3” per minute and, assuming observa- 
tions of seven-minutes’ duration, about ten observations can be made 
on a standard 18-exposure film. Each record is positively identified at 
the time recording is started by telegraphically coding the fiducial 
light source. 

To facilitate operation a hand hole is provided in the top of the 
pressure housing for removing the camera, changing lights, and mak- 
ing other routine adjustments. The hand-hole cover can be removed in 
less than one minute, thereby avoiding breaking of the main seal for 
all but major adjustments. 

A water alarm comprising a ball float in a sump is provided to 
inform the operator if water is getting into the container. This ball 
float closes a contact which turns ona red light on the right side of the 
control panel. The instrument is further protected against damage 
from leaks at both the hand hole and main seals by a spun aluminum 
water-deflector which fits inside the upper half of the pressure housing. 
Holes are drilled in the lower half of the housing to permit any leakage 
water to run to the sump in the bottom without dropping on any of 
the mechanism. 

The gravimeter used for the under-water unit is a standard Gulf 
instrument in almost every respect and can be converted to a land 
unit in less than one-half hour by merely taking it out of the pressure 
housing and placing it in a land-conversion mounting which includes 
a different top plate and a mounting assembly with hand-operated 
leveling screws. The main deviation from the standard land unit is 
the thermal insulation which is about two times as efficient. While 
this improvement provides a relatively small decrease in over-all bulk 
of a land instrument, a net saving of 200 to 300 lbs. was made in the 
weight of the underwater instrument by this change. 

The underwater assembly is constructed mostly of aluminum, 
weighs about 300 lbs. out of water and 25 or 30 lbs. under water. By 
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making the unit almost buoyant it is possible to place additional 
weight where it is needed to produce the greatest stability under 
water. Three hundred and fifty pounds of lead weights are fastened 
to the tubular supporting structure for under-water use. This addi- 
tional weight can be quickly removed when it is desirable to carry the 
instrument to a land subbase. It can be easily carried by three men 
when the weights are removed. 

The pressure housing is strong enough to withstand the collapsing 


(RETOUCHED) 





Fic. 4. Illustrating a record taken in quiet water. Note 
the telegraphic identification on the left. 


pressure of a 700 foot head of water, is made of an aluminum alloy, 
bolted together with stainless steel bolts and corrosion protected by 
zinc washers. 

The control unit contains seven meters, four pilot lights, nine 
switches of various descriptions, two relays, and two motor driven 
contactors for fine adjustment of the level mechanism. Meters are 
mounted horizontally in line on the top of the control unit to give the 
minimum disturbance with rolling and pitching of the boat and are 
covered with a sealed-in plate glass window to protect them from water 
damage. Switches and other controls are on the vertical front panel 
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4th. Reflection Retouched ) 





Fic. 5. This record was taken in rough water and is considerably poorer than 
the average record although it can still be read to 0.1 to 0.2 mg. 


which is protected by a cover sealed with a rubber gasket when not 
in use. 

There is positive indication in the control unit of every step and of 
the operability of every piece of apparatus involved in taking an ob- 
servation so that the operator has full knowledge of instrumental 
conditions at all times. Thus inadvertent damage to the apparatus 
due to the remote control operation is minimized. 

‘Power is supplied to the control panel by two 12-volt storage 
batteries. An additional 12-volt battery is used to maintain the tem- 
perature of the gravimeter, which is controlled to a thousandth of a 
degree over short time intervals and a few hundredths of a degree over 
long time intervals. 

In order to handle all of the controls and safety devices, a nineteen 
conductor electrical cable, V, Fig. 3, connects the control cabinet to 
the gravimeter. Four of the conductors are used for a common ground 
but provide one or more spare conductors in case of damage to other 
wires; the remaining 15 conductors control: 4 motors (3 of which are 
reversible), 3 separate light-source circuits, thermostatic heat control, 
15°-limit control, water alarm, clamped signal, and two photo-cell 
circuits which do not use the common ground. 
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Fics. 6 and 7 are pictures taken of actual operation in the Gulf coast. In Fig. 7 
the instrument is being lowered into shallow water in order that it may be carried to 
a land subbase to tie water data to adjacent land data. The hoist assembly operated by 
power take-off from the boat engine is clearly shown. 
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There are no relays in the gravimeter housing, all circuits being 
controlled directly from the control unit. There are, however, a total 
of three relays in the setup, one for the heat control, another for the 
clamped-unclamped indicator and the third for the 15°-limit control 
all in the control unit and associated parts. 

Records are read by a special micrometer slide-comparator in 
much the same fashion as the gravimeter is read directly by a mi- 
crometer eyepiece. The advantages of continuous photographic re- 
cording for this work are obviously its simplicity and the permanence 
of record, but more important is the fact that even though the instru- 
ment never comes to rest, a fairly good gravity value can be obtained 
by special reading procedure. Records are corrected for film shrinkage 
errors by reading the fiducial lines previously mentioned and correct- 
ing the values accordingly. 

Fig. 4 shows a typical record in quiet water. The heavy lines are 
fiducial lines. Line marked Z is the zero line. Relative gravity can be 
computed by measuring the distance between the “zero” and any line, 
R, and multiplying by the proper conversion factor or constant. This 
record can be read with an accuracy of about .o2 milligals. Whether 
or not the gravity value is good to that precision will depend to a 
great extent upon the drift curve. 

Fig. 5 illustrates a record taken in rough, shallow water. Proper 
reading of this record should give a gravity value with a probable 
error of 0.1 or 0.2 mg. 

A maximum of 17 stations per day has been taken to date with an 
average of 6 stations per working day. The principal causes of delay 
in this type of work are unfavorable weather conditions and the slow- 
ness of boat transportation. 

Another handicap due to the relatively slow speed of boats is the 
time consumed in travel between the scene of operations and head- 
quarters. While this handicap can be overcome to some extent by 
providing mobile quarter-boats for the crew, it is evident that work of 
this type can never attain the speed of land operations. 
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THE EFFECT OF DENSITY ON 
SEISMIC REFLECTIONS* 


S. S. WEST 


ABSTRACT 


By means of existing data on the density of rocks, it is shown that reflection coefii- 
cients of seismic waves, calculated from velocities alone, often differ greatly from the 
actual values. Data on a large number of density discontinuities have been tabulated, 
and the true reflection coefficients calculated for six cases in which velocities are also 
known. 


The reflection of sound waves at discontinuities in the earth, upon 
which the art of seismic prospecting is based, has nearly always been 
treated as though only the discontinuous change of the velocity were 
of any importance. The assumption was convenient in that it freed 
the geophysicist of the necessity of much tedious measurement, but 
it was certainly not a good approximation. The purpose of this paper 
is to point out that changes of density can be quite as important as 
changes of velocity in the production of seismic reflections and that 
therefore the empirical correlation of reflections with abrupt changes 
in the velocity-depth curve may sometimes lead to serious errors. 

The ratio R of the amplitude of the reflected wave to that of the 
incident wave for normal incidence at a plane interface is given by the 
well-known formula 


p2V2 — piVi 
p2Vo+ piNi 


where p is the density of the medium and V the velocity of the longi- 
tudinal wave. The index 1 refers to the medium in which the incident 
and reflected waves are propagated (upper stratum in our case) and 
the index 2 to the medium beyond the interface (lower stratum). The 
sign of this expression gives the phase of the reflected wave, for, when 
p2V2 is greater than piV;, a compression is reflected as a compression, 
and, when pV; is greater than p2V2, a compression is reflected as a 
rarefaction. Thus the first motion of the ground at the seismometer 
will be up when p2V2 is greater and down when p,V; is greater. The 


* Presented at the Annual Meeting, Chicago, IIl., April 9, 1940. 
} Joint Geophysical Laboratory of Stanolind Oil and Gas Company and Western 
Geophysical Company. 


45 








46 S. 8. WEST 


assumption of normal incidence is close enough to actual conditions 
to be a satisfactory approximation. 
The longitudinal velocity V is given by 


vw (- +t “\" 
p 


where A and yp are the Lamé constants. If we assume a fixed value of 
Poisson’s ratio o, then A is proportional to w and V=k(u/p)"/?, where 
k is a constant. This assumption is justified by the fact that o changes 
slowly in the earth, approaching the value j at large depths. In this 
latter case we may write the formula for the reflection coefficient R as 





= (poue)1/? sos (pipa)!/? 
(pope) !/? + (pimr)}/? 


where p is the rigidity of the medium in question. 

Although we have usually to do with p and V in practice, a classifi- 
cation of discontinuities with regard to p and m gives a clearer idea of 
the geological conditions which define them. In this system the fol- 
lowing eight types yield reflections. 

1). p2>p1, M2>pMi (V2 nearly the same as Vj). This is one type of 
true density discontinuity when V2=V1, ue/p2=mi/pi. It may occur 
when induration increases both density and rigidity. First arrival is 
upward, i.e., phase is unchanged by reflection. 

2.) po> pi, M2="i(V2<Vj). Probably not a very common case and 
not productive of large reflection coefficients when it does occur. 
First arrival is upward (phase unchanged). 

3+) p2> pr, Me<Mi(V2< Vi). This yields an unusually large velocity 
discontinuity, but a small reflection coefficient. It can occur with salt 
on anhydrite, or possibly in a homogeneous medium in which there is 
a zone of saturation by fluid which increases the density of the satu- 
rated medium. First arrival upward or downward, or no reflection. 

4.) p2=pi, M2 >Mi(V2>Vi). The commonly assumed true velocity 
discontinuity. First arrival upward (phase unchanged). 

5-) p2=pi, Me<pi(V2<Vji). The inverse of type 4. First arrival 
downward i.e., phase reversed by reflection. 

6.) po<p1, de>, (V2>Vi). The inverse of type 3. First arrival 
upward or downward, or no reflection. 

7.) pe<pi, M2=M1, (V2>V31). The inverse of type 2. First arrival 
downward (phase reversed). 
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Fic. 1. Portions of density logs which contain discontinuities of density sufficiently 





large to produce strong reflections or to modify greatly reflections from discontinuities 
of velocity at the same depths. 


I. 


qu 


8.) p2<pi, Mo<pi(V2 nearly the same as Vj). The inverse of type 
First arrival is downward (phase reversed). 

The data needed to illustrate these types of discontinuities are 
ite rare and hard to procure, but useful information can be obtained 


from several sources. Mr. R. C. Coffin of the Stanolind Oil and Gas 
Company has kindly provided density logs of ten wells located near the 
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TABLE 1 
Mayor Density DISCONTINUITIES FROM MEASUREMENTS ON CUTTINGS FROM WELLS 
Density 
Location of Well Depth R (%) 
Upper Lower 
Pure Oil Co. #1 Hendricks 2300 2505 | 2.78 12.8 
Winkler Co., Texas 2440 2.92 2.28 —12.3 
Texas-Rhodes 1140 2.87 2.12 —15.0 
Lea Co., N. M. 1220 2.07 2.78 TA: 7 
1600 235 | 2.10 —13.4 
1620 7 2.85 TAGS 
| 
1640 2587 2330 — 9.8 
1720 2:22 2.76 9.0 
1750 2.18 2.65 0:57 
1940 2.11 277 E3i5i5 
1980 2.18 2.66 9-9 
2000 2:62 | 2.13 —10.3 
2265 2.2) | 2.74 9.8 
2300 2.45 | 2.89 8.2 
2390 2.82 | 2.22 —= DEO 
2460 203 | 2.84 | 14.3 
2500 2.82 | B12 | 12.2 
2685 27 | 2.60 g.o 
Carter-Zweibel 1460 257 | 4:23 — 9.2 
Andrews Co., Texas 1565 2.23 2.98 14.4 
1870 2:43 2.03 9.3 
1895 2.93 2.45 — 8.9 
1985 2.18 2.63 9-4 
2195 2217 2.87 13.9 
2235 2 EF 2.94 1 ae | L 
2760 2.08 2.82 B51 
2780 2.82 223 —13.7 
Marland-Zock #1 Guitar 1215 2.80 2.28 — 9.8 
Howard Co., Texas 1700 2.33 2.93 II.4 
2270 2:20) 2.80 9.8 
Humble #1 Slaughter 2035 22 a7 12.4 
Martin Co., Texas 2060 2.76 2.13 —12.9 
2605 2.20) 275 F1:.6 
2980 2.32 2.76 8.7 
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TABLE 1 
Density 
Location of Well Depth R(%) 
Upper Lower 

Texas-Lockhart 1660 | 2.92 2.32 —II.4 
Lea Co., N. M. IgIo 2.96 2.42 —10.0 
1990 2.42 2.95 9-9 

2020 2.95 2.32 —12.0 

2410 2.68 2.42 — 9.4 

2620 2.38 2.82 8.5 

2650 | 2.94 2123 —12.9 

2780 | 2.32 2.74 | 8.3 





southeast corner of New Mexico. The density of the cuttings of these 
wells was carefully determined for each five-foot or each ten-foot sec- 
tion to depths of three to five thousand feet, the measurements being 
made by means of a Jolly balance with the sample submerged in oil. 
Seven of these logs show large variations in density, and from these 
a reflection coefficient R= (p2—1)/(p2+p1) was computed for each dis- 
continuity. (By assuming equal velocities to calculate this reflection 
coefficient we depart somewhat from physical reality, but the coeff- 
cient thus obtained is a convenient measure of the magnitude of the 
discontinuity of density. The true reflection coefficient is calculated 
later for those cases where the velocities are known.) In these seven 
wells there were 112 discontinuitities with | R| >5%, 63 positive and 
49 negative. The data on those with | R| >8% are listed in Table 1. 

In some of the cases listed, the thickness of the bed below the dis- 
continuity might not be great enough to permit a good reflection, but 


TABLE 2 
DENSITIES FROM GRAVIMETER MEASUREMENTS 




















Location Formation Density 

1. Caddo and Washita Quartermaster (red shale and sandstone) 2.4 
Counties, Okla. Cloud Chief (gypsum) 2.7 
Day Creek-Whitehorse (dolomite) ee 

2. Carter County, Okla. Trinity sandstone nas 
underlying Permian 2.4 

3. Llano Estacado above escarpment | 2.0 
escarpment, N. M. below escarpment | 2.4 
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in Fig. 1 are shown portions of the logs in regions where the formations 
are thick and the density differences large. The boundaries of the 
salt section are shown on the logs in Fig. 1B and 1C. Data from 
“‘well-shooting” in Gaines County, Texas, (about 50 miles away) indi- 
cate that the velocity in the layer above the salt is about 12,000 
ft./sec., in the salt about 15,000 ft./sec., and below the salt about 
11,000 ft./sec. 

From a paper by L. L. Nettleton! we take the densities of adjacent 


TABLE 3 
REFLECTION COEFFICIENTS CALCULATED FROM BOTH DENSITY AND VELOCITY 








| Densities| Velocities Ref ection Co- 
efficient 





ae | ua | Lybe of 
Location Formations |————— | —| Discon- 
Equal |"! | tinutt 

p2 Pl V2 Vi Actual Pir y 

















| 
| 

Texas-Rhodes | Anhydrite 2.12|2.87|15 ,000|12,000—3.9% 11.1%|rt.2| 6 
Lea Co., Salt | | | 
| 


N. M. 





Salt Anhy- 2.60|2.17|11 ,000|15 ,000| —6.4 —15.4 |o.6} 3 
drite | | | 











Pure Oil Co. #1 | Anhydrite 2.52/2.85|15,000112,000 EO EGoE |ns4| 16 


Hendricks Salt | | | 



































Winckler Co., |————————- i——| —- 
Texas Salt Anhy- |2.78/2.15 11000115 ,000| — 2.5 —15.4 |0.7| 3 
drite | | | 
Caddoand Quartermas-|2.7 2.4 9,000] 8,300] I0.0 AsO. |1.2 I 
Washita Cos., ter Cloud 
Okla. Chief 
Cloud Chief |2.1 2.7 11,000] 9,000)—2.5 4.0 |t.2| 6 
Day Creek | 














strata which are listed in Table 2. For the first three formations it is 
possible to estimate velocities from the shot point seismometer breaks 
of seismograph records made in that region. Where the formations 
are deeply buried the velocities will be greater, but their ratios (and 
the ratios of the densities) will probably be much the same. The veloc- 
ity thus obtained for the Quartermaster formation is 8300 ft./sec., 
that for the Cloud Chief is 9,000 ft./sec. and that for the Day Creek 
is 11,000 ft./sec. Of the other two pairs of formations, it can only be 


1“TDetermination of Density for the Reduction of Gravimeter Observations”; 
Geopuysics, Vol. IV, No. 3 (1939) p. 176. 

















EFFECT OF DENSITY ON SEISMIC REFLECTIONS 51 


said that the assumption of equal velocities gives them reflection 
coefficients of 7% and 9%, respectively. 

Six discontinuities are sufficiently well described by the above 
data to permit the calculation of the true reflection coefficients and 
comparison with reflection coefficients from the assumption of equal 
densities. Table 3 contains the data and the results of these calcula- 
tions. The ratio w2/u1, which is necessary for the classification of the 
discontinuity, is given by 


7,9 
J 2°p2 





M1 Vi"p1 

Thus, in five of the six cases we have been able to investigate the 
effect of the change of density has been to reduce the amplitude of 
the reflection so much as to make it scarcely detectable, and in the 
first and last cases the phase is reversed. One case (fifth line of Table 3) 
is a true density discontinuity of type 1, which type should be fairly 
common, since rocks often become more rigid as they become more 
dense, and conversely. The occurrence of types 1 and 8 especially needs 
investigation, for they produce reflections very difficult to correlate 
with velocity data. 

It will therefore be seen that one must take account of density as 
well as velocity in any quantitative consideration of seismic reflec- 
tions. The expense of obtaining density data should not be great where 
access can be had to cuttings from one or two wells in the neighbor- 
hood of a prospect, and geological and velocity data can be obtained 
from the same wells. In this way reflections can be assigned to geologi- 
cal interfaces much more definitely than has hitherto been possible. 
Further research on the problem should yield results of considerable 
value for the interpretation of the field data of seismic prospecting. 











RESOLUTION CONTROL IN SEISMIC SURVEYS* 
ROLAND F. BEERSf 


ABSTRACT 


Development of the concept of velocity stratification by experimental means has 
revealed the origin of certain reflection horizons. Correlation of these data with sub- 
surface lithology suggests a method of attack, employing seismic means, on the problem 
of the stratigraphic trap. 

This paper outlines the problem and describes a method whereby stratigraphic rela- 
tionships may be measured. It is suggested that a seismic wave may be fitted to the task 
by altering its wave length in accordance with the degree of resolution desired. Experi- 
mental examples of the procedure are given. 


Field experiments tending to develop the concept of velocity strati- 
fication have disclosed details of sedimentary stratigraphy which are 
of value in allocating sources of seismic reflections. As mentioned in a 
previous paper! these reflections are found to originate in stratified 
media where a sudden change in seismic properties is encountered in 
the path of a seismic wave. These discontinuities may be of one or 
more of several types. Included among them will be found sudden 
changes in density, modulus of elasticity, velocity of propagation to 
longitudinal waves, and acoustic impedance. These properties of 
rocks im situ are susceptible to direct measurement by experimental 
means. The results of these measurements yield stratigraphic logs 
which are valuable for the purpose of correlation between different 
locations where bore holes have been drilled and analyzed experi- 
mentally. 

In addition to the function of correlation these measurements yield 
detailed information regarding the nature and properties of rocks en- 
countered in the path of the seismic waves. This paper is primarily 
concerned with the distribution of these properties throughout the 
geologic column and the discussion of problems of correlation for 
which a solution may be available as a result of these studies. 

One of the most frequently occurring interfaces where good quality 
reflections have been observed is that of a limestone overlain by shale. 
The magnitude of seismic discontinuity at such a plane as determined 
from measurements will produce a computed reflection coefficient 


* Presented at the Annual Meeting, Chicago, IIl., April 9, 1940. 
} The Geotechnical Corporation, 1702 Tower Petroleum Bldg., Dallas, Texas. 
1 Roland F. Beers, GEopuysics, Vol. V, No. 1, (1940), pp. 15-21. 
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between 0.15 and 0.30; a reflection of this order is usually an outstand- 
ing or “dominant”’ reflection of a record. Reflection coefficients of lesser 
value down to o.10 are discernable and frequently offer support to the 
correlation pattern in their interval relationship to other “dominant” 
reflections. Supporting theoretical considerations, it is not infrequent 
that good quality reflections are found at the interface separating a 
shale overlain with limestone. In general, coefficients of this class have 
not been observed to be as large as those of the inverse sequence men- 
tioned above. There is an interesting relationship between two con- 
secutive reflections, the first from a shale to limestone interface, the 
second from limestone to shale. Examples of this reflection pattern 
are frequently observed in beds of Devonian to Ordovician age, such 
as are exemplified in and near the Seminole Plateau District of Okla- 
homa. This province is characterized by dominant reflections from or 
near the top of Hunton Lime, others within or near the base of the 
Hunton Lime if the limestone section is thick enough, and an addi- 
tional outstanding reflection from the interface between the Viola 
Limestone and the overlying Sylvan Shale. For many years mapping 
of structural details in this province has depended upon the ability 
of the seismic interpreter to follow accurately the correlations of re- 
flections from each of these stratigraphic positions. Not infrequently 
the intervals between adjacent reflections in this section undergo wide 
changes. These changes may have been caused by changes in the 
amount of deposition of the formations or by removal by erosion from 
the top of the formations. Intimately related to the problem is the 
influence of the Woodford or Chattanooga Chert which is frequently 
the beginning of a series of reflections of equal dominance spaced at 
equal intervals from one another. 

The approach to this intricate stratigraphic problem has often 
been by the route of fortuitous circumstances, although Morris Mus- 
kat has shown analytically the relationships which control the reflec- 
tion wave patterns in the three-layer problem.” It has been observed 
by many that the ability of the reflection seismograph to resolve the 
stratigraphy of this province into an orderly array has been affected 
profoundly by the circumstances surrounding the shot hole. In many 
instances reflection correlation profiles have been repeated again and 
again with the hope of obtaining just that particular character of rec- 
ord which would enable one to “see” clearly the exact stratigraphic 


2M. Muskat, Gropuysics, Vol. III, No. 3, (1938), pp. 198-218. 
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picture. One common practice employed in the field for influencing 
the situation has been the depth of shot hole. In some cases it has been 
found that a shallow hole produces more intelligible results than does 
a deep hole. On other occasions the reverse has been observed. Faust® 
has made observations on this point. 

A simple but valuable concept leading toward an understanding of 
this problem may be drawn from the field of engineering. Here it is 
common practice when making a physical measurement to select a 
measuring device consistent with the dimensions to be measured. In 
field surveys where a precision of a few feet is satisfactory this prac- 
tice leads to the use of the transit with stadia wires for dimensions of 
the order of a fraction of a mile. For finer resolution a steel tape is em- 
ployed whereby a distance of 100 feet may be measured within a few 
hundredths of a foot. For instruments of great precision micrometers 
are employed leading to an accuracy of .oor of an inch. In approaching 
the problem of stratigraphic analysis it is helpful to consider the 
length of measuring devices which we might employ. It will be seen 
that this choice falls naturally upon the wave length of the seismic 
wave which is to be fitted to the dimensions of the stratigraphic 
problem. 

It has often been customary to speak of seismic waves as having 
some definite frequency. Considerable effort has been expended in 
maintaining these waves as recorded on seismograph records within 
certain frequency limits. It seems likely to be more helpful for these 
purposes if we abandon this concept and view the reflection phenome- 
non in terms of the length of time during which the seismic impulse 
persists. This period of time multiplied by the instantaneous rate of 
wave propagation through the formation then determines the length 
of the wave which is at that instant available for measurement of 
stratigraphic intervals. If we wish to control dimensions of this seismic 
yardstick we have one variable upon which to operate: the duration of 
the seismic impulse. It seems possible to control the length of this 
impulse most effectively at its point of generation, i.e., the shot point. 
Certain it is that much control can be effected here with beneficial 
practical results. In the East Texas Basin, for example, the length of 
this impulse can be varied between limits of 0.02 and 0.04 seconds. This 
is common knowledge to a great many operators and is repeated here 
only in order to summarize pertinent information in one place. The 


3U.S. Patent No. 2,156,624, May 2, 1930. 
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manner in which this impulse can be varied is determined somewhat 
by the local conditions. These involve the nature of the soil, the depth 
to the water table, and the bed rock at the base of the aerated zone. 
Considerable effect is also imposed on the recorded seismic wave by 
the surface and aerated zone, but these effects are separate from 
those upon the generated wave and will be considered later. 

It has been observed in East Texas that a shot located below the 
top of the unweathered zone, also below the lower surface of ground 
water, will generate an impulse whose duration is close to the lower 
limit mentioned above. If the shot is pulled up into the main body of 
ground water the duration of the impulse will be lengthened, and the 
longest impulse of all will be generated when the shot is located just 
above the top of the water stratum, usually within the weathered 
zone. Unfortunately a shot located in this last mentioned position 
produces an undesirable ground wave of long duration and large ampli- 
tude. This wave, commonly spoken of as “‘ground roll,” travels hori- 
zontally along a path intervening between the shot point and re- 
corder stations and frequently succeeds in completely masking all 
reflection waves arriving coincident with it. It is customary in such 
circumstances to space the recorder and shot-point stations at such 
intervals that reflections may be observed on the record at times when 
no ground roll is present. In other instances the horizontally traveling 
wave is diminished in amplitude by the use of high-pass filters or 
through the use of interference effects between multiple detectors. 
None of these factors have any bearing on the problem under discus- 
sion but they are mentioned in order to eliminate effects which are 
often considered of importance in handling the present problem. It 
will be seen that they are irrelevant. 

The effect of altering the duration of the impulse generated at the 
shot hole will depend upon the stratigraphic problem under examina- 
tion. In East Texas one of the outstanding reflections commonly em- 
ployed for correlation purposes is the so-called “‘Georgetown” or 
“‘Comanchean.”’ Apparently this reflection comes from more than one 
interface in the geologic column. When the duration of the generated 
impulse is short, i.e., of the order of .o20 seconds, it is possible to ob- 
serve the component reflections, one coming from each surface of 
discontinuity. It will be seen that if the generated impulse could be 
adjusted to a certain critical length then the reflected components 
could be made to appear to be one single reflection of extended dura- 
tion. This process is accomplished by arranging to have consecutive 
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impulses of such duration and at such intervals as to combine addi- 
tively in phase, thereby producing a larger amplitude in combination 
than is possible either as individual pulses, or combined in any other 
phase relationship. This concept has been supported repeatedly by 
experimental means. The correct procedure in such a problem depends 
upon the aim of the survey. If it is desired to map structural details 
on the top of the Comanchean section without investigating the more 
minute intricacies of the stratigraphy then the most desirable impulse 
will be that of long duration. By its use each of the component reflec- 
tions will add together in favorable phase relationship so as to produce 
an outstanding reflection of larger amplitude than would be available 
if the components were separated. For this reason it is common to 
observe that shallow shot holes in East Texas often produce the clear- 
est results and that many prospects where deep holes have been drilled 
have been found to be incorrect in their structural interpretation. 

If impulses of short duration are employed, the confusion attend- 
ant upon correlation efforts becomes very great. This difficulty is 
often derisively spoken of as the “pulse” problem or the ‘‘leg”’ prob- 
lem. The number of “legs’’ with which one has to contend in making 
an interpretation under these conditions is so great that it is fre- 
quently impossible to stay on the same “leg’”’ of a reflection. As a 
result, mis-correlations frequently follow of the order of at least one 
“leg.”’ In East Texas this may amount to about 125 feet. 

In attacking the stratigraphic problem of the Devonian to Ordo- 
vician section of Oklahoma a different situation is often found. Here 
the ability of the seismograph to follow stratigraphic changes accu- 
rately calls for a fine degree of resolution. It is often necessary to know 
whether or not a reflection is being returned from the Hunton lime. 
In areas where thick Hunton is present it is necessary to know whether 
or not two Hunton reflections may be expected. Practical experience 
has shown that when the Hunton exceeds a thickness of 275 feet it is 
reasonable to expect one from the Hunton top and one from an inter- 
val of 275 or more feet below it. It will also be observed, following this 
sequence, that there is a reflection from the Viola, which, in the case 
of very thick Hunton, may be of diminished amplitude. This is be- 
lieved to occur because of the reduction in the amplitude of the seismic 
wave transmitted through the thick Hunton section to the Viola top. 
Frequently an estimate of the presence of the Viola reflection can be 
tentatively based on this consideration. Under these conditions it is 
desirable to employ a seismic impulse of as short duration as possible, 
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consistent with all other requirements. In general the shorter the 
duration of impulse employed the more intricate becomes the problem 
of correlation. So many ‘“‘legs’”’ can be developed on each of the reflec- 
tions with excessively short impulses that it becomes an intellectual 
impossibility for the interpreter to keep everything straight. There 
are so many possibilities of mis-correlation that he becomes hope- 
lessly confused. 

This is especially true if intervals occur which change; this is the 
problem of the stratigraphic trap. To reduce this difficulty, resort has 
been had in the past to continuous profiling. By this means correla- 
tions are established over subsurface intervals of fifty feet in horizontal 
extent. On the assumption that no extreme changes in depth to a re- 
flecting horizon can occur within this short distance, it is possible to 
make correlations on continuous profile records from trace to trace. 
Contiguous records are tied together by at least one overlapping posi- 
tion so that no difficulty is had in correlating from record to record. 
The only remaining difficulty is then to follow a reflection or series of 
reflections accurately across each record. If good quality records are 
obtained and if the resolution in use is well adapted to the strati- 
graphic problem, then accurate results may be expected. 

In this province a good shot hole is said to be one which yields 
sufficient but not unnecessary resolution. Fortunately this has been 
the prevailing situation here since the inception of the reflection seis- 
mograph. It is interesting to speculate upon the probable develop- 
ment of the seismograph if these fortuitous conditions had not pre- 
vailed in this province. This area was literally the cradle of the reflec- 
tion seismograph; under these ideal conditions many results were ob- 
tained of great practical value. This encouragement undoubtedly led 
to extension of the method to other less favorable areas where progress 
was made possible only through the faith engendered in proven pros- 
pects. 

There are many stratigraphic problems which do not yield solu- 
tions as easily as the one mentioned above. In some cases it is neces- 
sary to lengthen the duration of the seismic impulse generated at the 
shot point. In others it is desirable to shorten this impulse. Many 
workers have contributed to this problem and a number of patents 
have been granted on methods of control. Among those which have 
been examined by the writer are: 

U. S. Patent No. 2,151,878, B. B. Weatherby, 3/28/39, 
U. S. Patent No. 2,154,548, B. B. Weatherby, 4/18/39, 
U. S. Patent No. 2,156,198, S. A. Scherbatskoy, 4/25/39. 
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Other workers have undoubtedly observed this phenomenon and have 
contributed to its development. In general the principles employed 
are the same and are intended to produce a seismic impulse at the 
shot point of pre-determined duration. There is also some evidence of 
intention to shape the impulse, that is, to control the steepness of the 
wave front. This will be discussed in connection with receiving appara- 
tus and does not produce the same effect as does the control of actual 
pulse duration. 

Some have observed the influence of the dimensions of the shot 
hole cavity upon the characteristics of the generated wave. J. A. 
Sharp‘ has outlined the idealized conditions and has observed that 
the length of the generated wave depends, among other factors, upon 
the radius of the cavity created by the explosive. This radius is deter- 
mined by the physical characteristics of the shot hole medium. At the 
inception of the explosion a wave of spherical dimensions proceeds 
radially from its source. The explosive wave shatters the medium in 
its path so long as its force exceeds the “rupture strength” of the ma- 
terial encountered by the wave. The force diminishes. as the radius 
increases in a manner closely resembling the analogous radiation 
problem, due allowance being made for energy transformation into 
frictional and heat losses. The wave front finally arrives at some ra- 
dius where the medium is no longer shattered by the existing force, 
at which point propagation of the elastic wave takes place. It is the 
dimensions of the cavity at this instant which determine the wave 
length and other characteristics of the generated wave. These con- 
siderations are affected by water-loading of the shot hole medium and 
the effectiveness with which the shot is tamped. 

The dimensions of the cavity where elastic wave propagation sets 
in also affect the amplitude to be expected from reflected waves ob- 
served at the recorder stations. Assuming that the propagated wave 
may be considered ideally spherical, it is seen that the amplitude of 
the reflected wave for a given charge of explosive will be greater, the 
larger the radius of the cavity at the onset of the elastic wave propa- 
gation. 

In addition to the influence of the generated wave upon the control 
of seismic resolution, we have also to consider the characteristics of 
the seismograph recorder instruments. However these do not control 
resolution. They may prevent adequate control after it has been 


4 Unpublished manuscript, read at the S.E.G. Annual Meeting, Los Angeles, Calif., 
March 1937. 
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established. H. A. Wilson® has shown in his very valuable paper the 
response characteristics of seismographs whose natural period is 
longer than, equal to, or shorter than the period of ground motion. 
He has analyzed the response of these instruments when critically 
damped and when undamped. He makes the observation, ‘Critically 





Fic. 1. Response to two successive Heaviside unit-functions of 
opposite sign, separated by various intervals. 


damped seismographs are obviously better when it is desired to ob- 
serve the arrival of several waves separated by appreciable intervals.” 
This is clearly a contribution to the problem at hand and coupled 
with adequate control of the characteristics of the generated wave 
affords practical results which are otherwise unavailable. 

Fig. 1 shows a typical response curve taken from a critically 
damped instrument whose period was greater than the period of the 
exciting impulse.* The effect of lengthening the duration of this im- 
pulse may be noted from the latter parts of Fig. 1. As the natural 
period of the critically damped instrument is first equalled and then 
exceeded, it will be noted that the response recorded by the seismo- 
graph becomes increasingly complicated. These experimental results 


5 H. A. Wilson, Physics, Vol. 2, No. 3, (1932), pp. 186-199. 

* Note on the preparation of these oscillograms: A complete reflection seismograph 
was set up in the laboratory. The detector was buried in the earth. The moving system 
of the detector was given a displacement as shown on the first (top) trace of each 
figure, by means of an electromagnetic force. The magnitude and wave form of this 
applied force was controlled as shown in the Figures. The two lower traces show the 
response of a critically damped seismometer. 
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verify Wilson’s analytical derivations and emphasize the importance 
of fitting the period and damping of a seismograph wave detector to 
the problem in hand. 

These experimental results offer, incidentally, some insight as to 





Fic. 2. Response to two successive impulses of same sign. 





Fic. 3. Response to two successive impulses of opposite sign. 


the nature of the impulse which is commonly observed on the reflec- 
tion seismograph. From an examination of Fig. 2 it will be seen that 
a single unidirectional impulse produces a response not exactly like 
the reflection wave. When a positive upward deflection is combined 
with a negative deflection of equal magnitude, then we observe that 














RESOLUTION CONTROL IN SEISMIC SURVEYS 61 


the response of the seismograph approaches that observed in field 
practice. Records taken in these illustrations were made in the labora- 
tory under ideal conditions so that the results would be free from un- 
wanted interference from other waves. 
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Fics. 4 and 5. Response to two successive sine waves, separated by various intervals. 


Fig. 3 shows tracings of a more complicated impulse. This consists 
of a positive displacement of unit value followed in order by a negative 
displacement of equal value, then at the same interval by an additional 
negative displacement of unit value, and finally a positive displace- 
ment, returning the round to zero position. The wave recorded by the 
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seismograph for this type of excitation will be seen to be almost an 
exact reproduction of the customary reflection wave as observed in 
field practice. 

In order to illustrate the influence of control factors upon the de- 
gree of resolution attainable, these experiments have been performed 
with a critically damped instrument of fixed period ¢. The duration of 
the exciting impulse has been altered as well as the interval T between 
two consecutive impulses. Figs. 4 and 5 show the response of the seis- 
mograph to a wide assortment of impulse characteristics. An appre- 
ciation of the possibilities of the resolution control by these means 
will be best afforded from a close study of these records. 

As reported by Wilson, the critically damped instrument records 
quite accurately the second derivative of ground motion. Under some 
circumstances undamped instruments of long period perform simi- 
larly. If it is desired to record a series of consecutive impulses it is 
important that the period of the instrument and its damping be ad- 
justed in accordance with the duration of impulses to be received and 
the time interval separating them. It might be thought that an insight 
into this problem would be had more directly if the recorded reflection 
wave were simplified in its pattern. This would be true except for the 
fact that considerable advantage is gained through recording of the 
second derivative whose amplitude may be many times that of the 
record of actual ground motion. Enhancement of this property is 
often employed to advantage through steepening of the wave front 
of the generated wave which can be accomplished in many instances 
through a fortuitous choice of shot-hole medium. In these experiments 
the steepness of the wave front has been varied between limits of 
.005 and .oo2 seconds. Wilson has observed the importance of this 
effect and has predicted the value of a very short period, critically 
damped, seismograph. He was limited by the sensitivity of the me- 
chanical instrument, a fact which is now of no great importance on 
account of the amplification available through electrical means. 

One particularly undesirable arrangement is presented by an un- 
damped instrument whose period equals that of the ground motion. 
The response of such a seismograph to this type of excitation yields 
a series of almost continuous oscillations. 

Nearly as serious in its effect upon resolution is the critically 
damped seismograph whose period is shorter than that of the impulse. 
This choice consistently adds at least one extra “leg” on the end of the 
recorded wave, thereby adding to the confusion of the already over- 
worked interpreter. 
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In all of these experiments damping has been achieved by viscous 
means, using oil, or by electromagnetic means alone, whereby the 
degree could be varied over large limits, or by a combination of the 
two methods. The manner is unimportant. 

If the fundamental period and the damping of the reflection wave 
detector can profoundly affect the appearance of the recorded impulse 
then it is also true that the time constant and damping of the electrical 
circuits associated with the detector will have similar effect. In present 
day electric seismographs it is quite commonly the custom to discrimi- 
nate against unwanted waves by frequency selecting means. In some 
instruments preference for reflection impulses is achieved with the 
use of band-pass filters; in others through the use of high-pass and/or 
low-pass circuits. Due consideration must be given to the design of 
filter circuits as well as of associated networks which have potentially 
resonant characteristics. It is possible to achieve considerable altera- 
tion in the shape of recorded impulses by the use of networks of this 
class. In general their effect is quite similar to that pointed out by Wil- 
son and qualitatively at least his relationships may serve as a guide. 
These predictions have been verified experimentally by employing 
electrical networks having a wide range in values of damping and of 
tin.e constant. The results of these experiments are identical in char- 
acter to those illustrated for the earth detector shown in Figs. 3 to 5 
inclusive. 

It has been observed that the surface soil and aerated zone fre- 
quently exert some effect on the shape of the reflection wave. Many 
have pointed out from field work that if the detectors are placed in 
loose sand, records obtained from the same reflecting horizon cannot 
be correlated with those obtained when the detectors are placed in 
moist clay. Studies have been made intending to determine the influ- 
ence of surface materials and have been reported by Howell, Keen 
and Thompson.® These experiments were conducted with continuous 
waves and pulses and the writers report that insufficient work was 
done to achieve results of immediate practical value. Consideration 
must be given to the effect of this portion of the seismic reflection 
wave path in establishing control of resolution. Experiments similar 
to those reported by the writer for stratigraphic analysis! are under 
way and it is hoped that they may be the basis of a future paper on 
this important subject. 


6 GreopHysics, Vol. V, No. 1, (1940), pp. I-14. 








A NUMERICAL METHOD OF COMPUTING DIP DATA 
USING WELL-VELOCITY INFORMATION* 


KENNETH H. WATERSt 


ABSTRACT 


A method of computation for calculating steeply dipping strata is presented. It 
takes full advantage of the velocity information afforded by well shooting and is simple 
in application. Charts are derived from the velocity data which will give depth and posi- 
tion of the reflecting point together with angle of dip and azimuth of dip immediately 
from the time and differential time of a reflection taken from a seismic record. 


The ability to compute the angle of dip and the position of the 
reflecting point for a curved-path assumption has long been acknowl- 
edged to be almost essential in some regions, more especially in areas 

-where steeply dipping formations are encountered. Numerous at- 
tempts have been made and much ingenuity has been displayed in 
fitting algebraic laws to the change of velocity with depth, or some- 
times with time, and, for this initial assumption, in finding the path 
travelled by a ray. The type of function is limited however since any 
but the mathematically simplest cannot be integrated at a later stage 
in the analysis. Unfortunately these simple laws do not fit the well- 
shooting data very accurately. Mott-Smith' has shown that the rela- 
tions in which the velocity is expressed as a function of the time are 
more convenient to handle and fit the observed values better than rela- 
tions giving the velocity as a function of the depth. He also states 
that more complicated velocity-depth curves “‘such as the S-shaped 
ones” have not been investigated. It is the purpose of the present 
paper to show how all velocity-depth curves may be handled. The only 
restrictions which will have to be placed on the velocity variations 
are those which must be fulfilled in order that the variations are physi- 
cally possible, that is, that there must be a finite number of turning 
points and no infinities in the region to be examined. The velocity 
laws which will be studied have one common characteristic. The 
velocity is the same at any point in a horizontal plane at a given depth. 
This restriction will, of course, eliminate the computation of dips when 
the equi-velocity surfaces are anything other than horizontal planes. 


* Manuscript received May, 1940. 
+ National Geophysical Co. 
1M. Mott-Smith, Gropnysics, Vol. IV No. 1, (1939). p. 8. 
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It admits, therefore, of no change of velocity at a boundary between 
two different layers unless this boundary be itself horizontal. If such 
a boundary is near the surface and is known to be inclined, then its 
effect may be counteracted by applying a special form of ‘‘weathering”’ 
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Fic. 1. Diagram showing relationship between shot point and reflecting point 
and seismometer spreads on the ground surface. 


correction. If it is not near the surface, then other measures will be 
necessary. 

Except from an academic viewpoint, it is not necessary to know the 
path of the ray. All that is necessary in practice is a knowledge of the 
angle of dip of the horizon from which the reflection comes, and the 
position of the reflecting point relative to the shot-point. This position 
is usually given in cylindrical coordinates, the vertical through the 
shot-point being the cylindrical axis. The azimuth is referred usually 
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to the north or to some other convenient direction. In order to deter- 
mine the azimuth two seismometer spreads at right angles are taken 
(see Fig. 1). Then, since cylindrical symmetry is assumed for the 
velocity function 

ATr 


ATy 





tan ¢ = 


where ATg and ATy are the differential times of arrival at opposite 
ends of the east and north seismometer spreads respectively and ¢ is 


tin V-dV 








Fic. 2 


the azimuth of the dip. It is assumed, of course, that weathering and 
elevation corrections have already been applied to the differential 
times. 

It now remains to find the vertical depth (Z) and the horizontal 
displacement (R) of the reflecting point from the shot-point. It is pro- 
posed to adopt a numerical method for finding these coordinates. This 
numerical method will have no disadvantage from the strictly analyti- 
cal methods in accuracy since both of these rest finally on the accuracy 
of the original ‘‘well shooting” data. The mechanical integrations 
performed may be made as accurate as necessary by taking a suffi- 
ciently large scale for the graphs involved. It will be assumed that the 
seismic results comprise the total time (7) for the reflection and two 
differential times A7z and AT y from the two seismic spreads at right 
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angles. For simplicity, also it will be assumed that these spreads are 
distributed symmetrically about the shot-point and that one of these 
spreads lies in a north-south direction. 

Some formulae derived in this analysis have been obtained pre- 
viously by M. M. Slotnick? and have been quoted several times by 
other authors.’ Since, however, the method of derivation is illustra- 
tive of the mechanism of this method it will be given briefly below: 
Consider a thin layer (see Fig. 2) of thickness dz at a depth z below 
the shot-point. Let the seismic velocity of the layer be V and the 
velocity immediately above it be V—dV. Let a ray strike the upper 
surface of the layer, making 6—dé@ with the vertical. Then, by Snell’s 
law, if refraction takes place 


> 
‘ 


sin 6 | V 


sin (6 — dé) V—dV 





or 
cos 6d@ = dV. 


If this elementary equation be integrated it may be found that 


sin 6 Vz 
“eee (1) 


= 
sin 0 Vo 
where 6 is the angle to the vertical at a depth z and @ is the angle at 
the surface, V» the velocity at the surface, V, is the velocity at any 


depth z. 
In particular it should be noted that 


sin a Vr i 
si | (2) 


sin 0) : Vo 





where a@ is the angle of dip, Vr is the velocity at the reflecting point. 
It may be proved by purely geometrical considerations that the differ- 
ential times from records taken on two spreads at right angles may 
be compounded to give the corresponding differential time along with 
the direction of maximum dip according to the formula 


AT = [(ATx)? + (ATw)?]!”. (3) 
2M. M. Slotnick, GEopuysics, Vol. 1 No. 1, (1936), p. 9. 


3 C. E. Houston, Geopuysics, Vol. 4 No. 4, (1939), p. 242. 
4C. H. Drx, Geopnysics, Vol. 4, No. 1 (1939), p. 81. 
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If X is the effective length of the seismometer spread used it is evident 
that 


V AT : 
= sin Oo, (4) 





and hence, utilizing the relation (1) above 


V AT 
= sin 0. (5) 





The total time for the energy to travel from the shot-point to the 
reflecting point and back to the seismometer spread is given by 


zds 
T = 2 ——S (6) 


o Vz 


where ds is an element of path. 


: €& ‘ dz 
T = of arent A of - a 
0 Vzcosé@ o V(r — sin? 6)!/? 


{ dz ' (7) 
2 . 
0 x AT y 


ie J ( ) 
r ie 


introducing the symbol Vy to represent the phase velocity along the 
ground surface (=X/AT7) this formula may be written 


z dz z 
T= of = af F dz. (8) 
0 I ui ats 0 
Vi mR = cos 
V2 Vx? 


In this way the important relationship is brought out that no reflec- 
tion will be received if its phase velocity along the ground surface is 
less than the velocity at the reflecting bed. For convenience the inte- 
grand of this integral will be denoted by F,. 

In a very similar way -the horizontal displacement (R) may be 


obtained. 
R -{ tan 6dz 
0 
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Fic. 3. Interval velocity curve on which sample work is based. 
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Fic. 4. Curves showing the relationship between F; and ; for various values of 
AT/X. Values of AT/X marked on the curves are in units of 107%. 
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4 sin 6 
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The integrand of this integral will be denoted by F2. Suppose now that 
AT is known from the records, the function F; can be calculated for 
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Fic. 5. Curves showing relationship between T and z. AT/X values 
on curves are given in units of 10°. 
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any depth since the interval velocity (V,) is obtainable immediately 
from Fig. 3. A set of curves may thus be drawn showing the relation 
between z and F,. A sample set of such curves for various values of 
AT/X are given in Fig. 4. 

The area between any curve in Fig. 4 from the z axis to the curve 
from the line z=o to any depth will then give half the total time for 
the wave path to a reflecting point having that depth and the angle of 
dip specified. Thus for any given differential time a curve may be 
drawn snowing the relation between T and z. The time, however, is 
known directly from the record as is also the differential time and this 
latter curve will therefore give a direct method of finding Z. Sample 
curves, calculated from the interval velocity curve (Fig. 3) are given 
in Fig. 5. The procedure in finding the depth (Z) is simple. The differ- 
ential time appearing on the record is taken, AT/X calculated, the 
corresponding curve in Fig. 5 selected, and from the value of T on the 
record, the corresponding depth is obtained. 

By plotting a set of curves showing the relation between F2 and 
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Fic. 6. Set of curves for various values of AT/X showing relation between 
F, and z. AT/X values are in units of 107. 
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z and integrating these mechanically in precisely the same way that 
the integration of the F; curves were effected, the relations between 
R and gz for various values of the differential times are obtained (see 
Figs. 6 and 7). The value of Z has already been obtained from Fig. 5 
and by selecting the proper curve (with the correct value of AT/X) 
the value of R is obtained. 
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Fic. 7. Curves showing variation of R with z for various values of 
AT/X. AT/X is marked on curves in units of 10-*. 


The last stage in the calculation is the determination of the angle 
of dip. Fig. 3 shows the variation of the interval velocity with the 
depth as determined by “‘well shooting.” Since Z is already known it is 
a simple enough matter to determine V, from this curve and thence 
to obtain the angle of dip from the formula 

AT 
sina = Vr:-—- (10) 
x 

Although the method sketched above will doubtless leave very 
much to be desired from the point of view of mathematical complete- 
ness, especially with regard to the actual path taken by the seismic 
waves, it nevertheless presents one solution of the problem of utilizing 
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the knowledge obtained from shooting a well in the determination of 
dips when the strata are dipping steeply. Hitherto such knowledge 
could be applied only if, without much error, the beds were dipping 
so slightly that they could be assumed flat. The advantage of the 
method proposed is that any velocity-depth law may be utilized no 
matter how far it may depart from simple mathematical formulation. 
The method is simple and, once the necessary curves are obtained, is 
quick in application. The most obvious disadvantage is the amount of 
labor involved in computing the values used in drawing the curves. 
This, however, is largely offset by the fact that it is necessary to make 
these laborious calculations only once for each well shooting, and the 
curves obtained may be used in hundreds of dip determinations. 








A TRANSFORMED WAVE-FRONT CHART* 
RAYMOND A. PETERSONT 


ABSTRACT 


In seismograph computations assuming a uniform rate of velocity increase with 
depth, a “‘wave front-wave path” chart consisting of sets of orthogonal circles is often 
used. The present paper describes a transformation by which the wave fronts become a 
set of concentric circles with the wave paths as radii, while the cartesian coordinates 
representing vertical and horizontal distances in the ground become orthogonal circles. 


As an initial approximation in reflection seismograph computa- 
tions, the velocity of wave propagation in the ground is often assumed 
to increase at a uniform rate with depth. That is, if v denotes the 
velocity at depth y below the shot origin, then: 


v= % + ay. (1) 


The equations of the wave fronts and wave paths for this case have 
been derived in papers by Ewing and Leet! and by M. M. Slotnick.?* 
A wave generated at the shotpoint travels outwards in the form of an 
expanding sphere (Fig. 1) which drifts downwards as the radius in- 
creases. At any time ¢ after the shot instant, the radius of the wave 
front is given by: 


R = K(sinh at), (2) 
where 
K = v/a, 
and the depth of the center below the shot origin by: 
H = K(cosh at — 1). (3) 


Likewise, the wave paths are circles passing through the shot origin, 
with their centers on the plane: 


* Presented at the Annual Meeting, Chicago, IIl., April 9, 1940. 

t United Geophysical Co., Pasadena, California. 

1. Maurice Ewing and L. Don Leet, “Seismic Propagation Paths.” Transactions of 
the A.I.M.E., Geophysical Prospecting, 1932. pp. 245-262. 

2 M. M. Slotnick, ‘On Seismic Computations, with Applications, I.’”” GEoPpHysIcs, 
Vol. 1, No. 1. (1936). pp. 9~22. 

3 M. M. Slotnick, “‘On Seismic Computations, with Applications, IT.” GEopHysiIcs, 
Vol. 1, No. 3. (1936). pp. 299-305. 
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a 
and with radii: 
Q = K csc %, (4) 
where 
6) = angle of emergence of the wave path, Fig. 1. 


The angle @ between the wave path and the vertical at travel time ¢ is 
given by the relation: 


tan 6/2 = (tan 09/2) -e'. (5) 


If a wave front impinges normally at point P (Fig. 1) on a reflect- 
ing interface dipping with angle 0, the path of the reflected wave will 
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Fic. 1 


coincide with the path of the incident wave. If the arrival time, 7,, 
of the reflected wave is measured on a surface geophone-spread sym- 
metrical about the shotpoint and oriented in the direction of dip of 
the reflecting interface, the angle of emergence, 9, can be calculated 
from the relation: 


sin 09 = vo(dT,/dx) (6) 


where 


x = the horizontal coordinate along the spread. 
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The dip, 8, can then be calculated from equation (5). 

For computing the position of the reflecting interface it is conveni- 
ent to construct a “‘wave front-wave path” chart showing the position 
of the wave front for progressive intervals of time, and of the wave 
paths for various values of 6) (Fig. 2). In the present case the chart 
will appear as two orthogonal sets of circles superimposed on cartesian 
coordinates representing vertical and horizontal distances in the 
ground. If the wave front circles are graduated in terms of reflection 
time at the shotpoint, and the wave paths in terms of sin 4, the posi- 
tion of the reflecting interface can be read directly by finding the 
intersection of the proper coordinates on the chart. 


¢ 
Zz 
A 
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WAVE FRONT - WAVE PATH 
CHART 


Fic. 2. Wave front-wave path chart. 


Fig. 2 shows schematically only a few of the lines necessary for 
accurate interpolation. In an actual chart with four superimposed sets 
of closely spaced coordinates, repeated reading of data is a tiring proc- 
ess. The purpose of this paper is to describe a transformation by which 
the wave fronts become circles concentric about the shotpoint, and 
the wave paths become radii through the origin. The wave fronts and 
wave paths can then be replaced on the chart by a mechanical arm 
pivoted at the shot origin, and graduated with a scale showing the 
wave front position for progressive intervals of time (Fig. 3). The sine 
of the angle of emergence, sin 6), may be read directly on a circular 
scale at the lower end of the pivoted arm. In the process of the trans- 
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The equation of the wave front circles may be written: 





formation, the cartesian coordinates are warped into orthogonal sets 
of circles (Fig. 3). The transformed chart is used in the same fashion 





as the original chart, except that the reflection time and sin 49 are 


read on the appropriate scales. 
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(« — Q cos 0)? + (y + K)? = Q%, 


Fic. 3. Transformed wave front-wave path chart. 


x? + [y — K(cosh at — 1) |? = K? sinh? af. 





The equation of the wave path circles may be written: 


or 
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or 
(x — K cot 0)? + (y+ K)? = K? csc? >. (8) 


If the substitutions: 
sinh at = csch U, 


cosh at = coth U, (9) 
M = V, 
are made, equations (7) and (8) become: 
x? + [y — K(coth U — 1)|? = K? csch? U, (10) 
(x — K cot V)? + (y + K)? = K? csc? V. (11) 


Now equations (10) and (11) are identical in form with those for 
the equipotential surfaces and lines of force in the electrostatic field 
between two charged, parallel, circular cylinders.4 Accordingly, the 
methods of conformal transformation commonly used in potential 
theory can be applied to the wave front chart. 

Let the original chart be represented in the z plane, where: 


s=a+iy 


/ 


ol 


(12) 





j 
and the transformed chart be represented in the ¢ plane, where: 
f= E+ fa. (13) 


It can readily be shown‘ that the quantities U and V defined by equa- 
tions (10) and (11) may be represented by the real and imaginary 
parts of the function: 


— 
— 
1 


U+jV 
In (2 + 2j/K)/z. (14) 


To accomplish the desired transformation of the wave front chart, let 
jK 
W = |ln—.- (15) 
f 
Separating (15) into its real and imaginary parts, and eliminating U 
and V, there results: 
‘ William R. Smythe, “Static and Dynamic Electricity.””, McGraw-Hill, 1939. 
PP. 74-76. 
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£2 + 9? = K?[(cosh at — 1)/sinh at]? (16) 
and 
t/n = tan Oo. (17) 


Equation (16) demonstrates that the transformed wave fronts are 
circles concentric about the shotpoint, (Fig. 4), with radii equal to: 





K(cosh at — 1)/sinh at = K(H/R). (18) 
SHOT POINT F 
y 
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Fic. 4. Transformed wave front chart. 


The wave paths as indicated by equation (17) are radii through the 
shotpoint, making an angle 0) with the vertical. 

The equations of the transformed cartesian grid in the ¢ plane 
may be obtained by equating (14) and (15). Then: 


§ = jK2z/(z + 2jK) (19) 
or 
3 = 2jK¢/(jK — $). (20) 


Separating the real and imaginary terms in (20) and regrouping terms, 
there results: 
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. K? + Ky\? Kk? 2 
i (21) 
2K + y 2K + y 


(K2/x)2. (22) 





and 


G~ djs) + iq — £P 


From (21) it follows that the transformed y coordinates are circles 
having their centers on a vertical line through the shot origin, and a 
common point of tangency at a distance K below the origin. The radii 
of the circles are equal to 


K/(2 + y/K). 


Likewise, from (22) it follows that the transformed x coordinates are ° 
circles having their centers on a horizontal line at a distance K below 
the shot origin, and a common point of tangency vertically below the 
origin. The radii of these circles are: 


K*/s. 
From the above data, a transformed wave front chart can be con- 


structed for any particular values of vo and a. If desired, the chart 
may be calibrated in terms of dimensionless parameters.* 


* M. M. Slotnick has pointed out that the transformation described in this article 
can be accomplished by the method of inversion. 











DISCUSSIONS AND COMMUNICATIONS 
TIDAL GRAVITY OBSERVATIONS 


Last August. I sent you the manuscript of a paper on tidal gravity variations. Since 
that time, we made another test run of the gravimeter with very similar results. I am 
sending you a chart of the curves of gravity variation which we obtained. 

Curve A in the chart shows instrument readings converted into milligals, at twenty- 
minute intervals. 

Curve B is the gravimeter drift curve, obtained by plotting twenty-four hour run- 
ning averages of instrument readings. 

The observed tidal gravity variation is shown by curve C, obtained by subtracting 
the drift from observed instrument readings. 

The broken line, curve D, represents the theoretical tidal gravity variation for an 
unyielding Earth. 
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Fic. 1a. Tidal gravity variation 
A—Gravity meter readings 
B—Drift curve 
C—Observed gravity variation 
D—Theoretical gravity variation 


I computed the value of the factor 1+/—3k/2 by comparing the maximum varia- 
tion of gravity for successive cycles of the observed and of the theoretical curves, C and 
D. The average value obtained in this way is 1.2, the individual factors range from 1.1 
to 1.4. This average value is the same as that obtained by a similar procedure from 
observations made in December, 1939. 
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Fic. 1c. Tidal gravity variation. 
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An interesting fact was brought out in regard to the drift of the gravimeter. The 
first part of the observations, up to September 5, 20* G.C.T., was taken without 
clamping the gravimeter between the readings. Between that time—September 5, 20° 
—and the end of the run the gravimeter was clamped after a reading was taken and un- 
clamped just before the next reading. An examination of curve B shows that the drift 
of the instrument was very small, less than the probable error of a reading, as long as 
the instrument was left unclamped. The drift of the gravimeter in the second part of 
the record—from September 5, 20" G.C.T. up to the end of the record—when the in- 
strument was clamped between readings, was very much greater and of the same order 
of magnitude as during the run made in December, 1939. 

October 4, 1940 
A. WOLF 
Geophysical Research Corp. 
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REVIEWS 


Exploration Geophysics, by J. J. Jakosky. Times Mirror Press, Los Angeles, Cal., 1940. 
786 pp. $8.00 


The author of this book has a substantial background of practical as well as some 
teaching experience in his subject. He is well known as a contributor to various geo- 
_ physical subjects, but especially to the art of electrical prospecting and his treatment 
of this branch of exploration is exceptionally well done. The chapters devoted to it are 
outstandingly the best in the book. The chapter on Gravitational Methods is the 
weakest. 

The entire book has a distinctly practical flavor. Jakosky is evidently interested 
in theories and principles primarily because they have actual or potential application 
in the business of prospecting. 

The first two chapters are written primarily for those who deal with geophysical 
methods administratively, it being assumed that they may do some reading on the 
subject. These chapters deal with the nature and application of various methods to 
petroleum geology and mining, and to various engineering problems, including those of 
water supply and storage, civil and highway engineering. 

There follow five chapters which deal with the four geophysical methods, the use 
of which has led to important discoveries of minerals. Two of these chapters deal with 
electrical methods and the others with the magnetic, gravitational and seismic tech- 
niques. These chapters present the underlying physical principles, the properties of the 
rocks which are involved and the various techniques which are available for applica- 
tion to the practical solution of geophysical problems. The instruments which have 
been developed and the problems of actual field operations are given due consideration. 
Results of actual surveys are presented and their interpretation is discussed. 

Literature references are given as footnotes. They are hardly sufficient to provide 
ready access to much of the pertinent available literature. A useful feature is the list of 
United States, British, Canadian, French, German and Russian patents given at the 
end of each chapter. The value of these lists is somewhat impaired by the fact that some 
patent references are given as footnotes in the body of the text and these apparently 
have not been included in the list at the end of the chapter. The repetition of these 
patent references to make the lists complete would have greatly enhanced their value. 

The remaining chapters of the book deal with geo-chemical and geo-thermal meth- 
ods, drill hole investigations, physical principles applied to production problems and 
with permit and trespass practice and insurance. The presence of the chapter on pro- 
duction problems in a treatise on exploration is somewhat surprising, but the echo 
method for determining the depth to the fluid level in wells is essentially a geophysical 
technique and the application has been developed by geophysical operators, Jakosky 
being one of them. 

The book reads easily because the author has the facility for presenting ideas 
clearly. There are occasional exceptions to this as for example, the description of the 
Hoyt gravimeter (p. 174). Jakosky’s Exploration Geophysics offers a satisfactory sum- 
mary of this important art brought down to date. It will be of interest to readers pro- 
fessionally engaged in this art, and it should be well received by the profession. 

E. A. ECKHARDT 
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Geophysical Exploration, by C. A. Heiland. Prentice-Hall, Inc., New York, 1940. 1013 
pp. $10.00. 


Professor Heiland enjoys the distinction of heading the oldest department of ap- 
plied geophysics in an educational institution in the United States. He has also dealt 
with the subject as a consultant and as a geophysical contractor. The author’s preoccu- 
pation with geophysics as an educational subject is strongly reflected in the pages of 
Geophysical Exploration. The book is encyclopaedic and academic in character. This 
is by no means to say that the practical aspects of the subject are neglected. The 
subject is comprehensively treated, in some cases exhaustively so. The discussion of the 
Eétvés torsion balance and directly related material is very complete and excellent, 
but the amount of space devoted to an instrument and related techniques which are 
practically obsolete violates the reviewer’s sense of proportion. The treatment of mag- 
netic instruments available for prospecting is excellent. 

The book is divided into two parts. The first part of 62 pages explains the nature 
and uses of the methods covered in greater detail in the second part of the book. In 
this part there are no equations. It is intended for readers who may be interested in 
knowing what the geophysical methods are and how and why they may be used. Brief 
chapters discuss the application of these methods to oil exploration, mining and en- 
gineering. In part two the four major methods of geophysical prospecting, gravita- 
tional, magnetic, seismic and electric are discussed in detail in as many chapters. At 
the beginning of each chapter one finds a definition of the symbols used, followed by a 
suitable discussion of pertinent physical theory. The physical properties of the rocks 
which are involved are discussed and laboratory and field methods for their determina- 
tion. Excellent tables give values for many rocks of interest to the practical geophys- 
icist. The special techniques of the method are discussed and the instruments available 
for their application are reviewed, In most cases illustrative examples of actual sur- 
veys are given, but in the chapter on Seismic Methods, this illustrative material is very 
meager. 

The two final chapters cover Geophysical Well Testing and Miscellaneous Geo- 
physical Methods. 

The reviewer is compelled to challenge the statement on page 19 that magnetic 
anomalies of over 50 gammas may be assumed to arise from variation in elevation or 
composition of the underlying igneous rocks and that anomalies of less magnitude are 
usually due to similar variations in the sedimentary rocks. On page 425 it is stated that 
the structure at Oklahoma City corresponds to a negative magnetic anomaly. Measure- 
ments made by reviewer’s company in 1928 after only the discovery well had been 
drilled show an orthodox magnetic high. 

Literature references are given as footnotes throughout the text and are reasonably 
adequate. The absence of any reference to a number of important papers which readily 
come to mind detracts from the value of the book in affording easy access to the perti- 
nent literature in any field. Reference to patents seems to have been studiously avoided 
and the few that have been mentioned are found in the body of the text where they 
would be hard to find if such references are sought. 

Dr. Heiland’s book is easily the most comprehensive one on the subject of explora- 
tion geophysics now available. It should find favor with all who are actively interested 


in its subject. 
I. A. ECKHARDT 
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Geophysical Prospecting for Oil, by L. L. Nettleton. McGraw-Hill, 1940. 650 pp., $5.00. 


In his preface, Dr. Nettleton states that his book is an outgrowth of a course of 
lectures given for several years at the University of Pittsburgh, and that the work is 
intended for the student and lay reader rather than for the geophysical specialist. The 
volume is nevertheless of considerable interest to the practicing geophysicist and will 
form a worth-while addition to his library. There has been a real need for a concise but 
comprehensive treatment of the principles of the various geophysical methods and their 
application to prospecting for oil. By using a direct approach to his material and avoid- 
ing verbiage, Dr. Nettleton has succeeded in meeting this need in a very satisfactory 
manner. 

The book is divided into four sections. One of 146 pages, deals with Gravitational 
Methods, another, of 74 pages, with Magnetic Methods. The Seismic Methods are dis- 
cussed in a section of 132 pages and the final section of 156 pages discusses Electrical 
and Miscellaneous Methods, together with measurements of various types commonly 
made in drill holes. A brief and sensible discussion of the interpretation of geophysical 
data forms a concluding chapter. The presentation of material in each section follows 
a logical system, beginning in each case with a discussion of the fundamental principles 
involved, and the relationship of the quantities measured to geological structure. The 
instruments used are then briefly described, and the field techniques employed are dis- 
cussed in some detail. The method of reducing observations is then outlined and usable 
systems of computation pertinent to the method are presented. Finally, the interpre- 
tation of the results obtained is discussed. Appended to each section is a well selected 
bibliography which, while not exhaustive, is a practical one for reference purposes. 

The section on Gravitational Methods is particularly good. The theory of the 
Eétvés balance is given and the precision pendulum is briefly discussed, but the author 
has very properly emphasized the present day importance of the gravimeter. The dis- 
cussion of this instrument and its use leaves little to be desired even by the practicing 
geophysicist. 

The section on Magnetic Methods points out the similarity of all potential field 
methods and so ties in very nicely with the preceding discussion of gravitational meth- 
ods. At the same time the differences between magnetic and gravitational effects are 
pointed out to explain the limited usefulness of magnetic methods. 

The section on Seismic Methods is well adapted to the needs of the student, but is 
of less interest to the geophysicist practicing in that field. The professional geophysi- 
cist may object to some of the author’s statements, even though recognizing that they 
are made in the interest of simplification. For example, the discussion of weathering 
corrections, and the weathering computation schemes presented in connection there- 
with, leave the reader with the rather definite impression that the so-called weathered 
layer is, in fact, a discrete surface layer of low velocity. Similarly, in a very brief men- 
tion of the application of curved-path wave theory to reflection shooting, the statement 
is made that the assumption of linear paths will rarely result in errors in the computa- 
tion of depth to a reflecting layer; but no mention is made of the possible errors which 
may be introduced into the computations of steep dips. On the other hand, the treat- 
ment of the subject is quite complete in that it outlines usable field techniques and defi- 
nite systems of computation. After perusing the chapter the student will have a very 
definite idea of “‘how it is done” although he can hardly expect to be able to do it. 

The final section of the book is rather sketchy because of the large number of 
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topics covered. Electrical prospecting methods are treated briefly but in sufficient de- 
tail to explain their relative unimportance in the search for oil. A chapter is devoted to 
a discussion of resistivity and natural potential surveys in oil wells. The failure to in- 
clude here a reproduction of a typical electrical log is rather startling since the book, in 
general, is extremely well illustrated. Geochemical methods are mentioned only in 
passing as “‘not proven.” Mention is made also of various methods used in wells such as 
core orientation, determination of fluid levels, and radioactivity logs. 

The author’s final chapter on Geophysical Interpretation will be interesting read- 
ing for all geophysicists since it includes a very sane view of the problems arising when 
physics meets geology. 

Dr. Nettleton’s book should be an excellent text for use either in a course designed 
for undergraduate physicists who intend to specialize in geophysics or in a course de- 
signed for geologists and petroleum engineers who wish to obtain an understanding of 
the various geophysical methods. As a textbook it has the very definite advantage of 
having been written by a practicing geophysicist, and from the physical point of view. 
The author has not hesitated to use mathematical treatment wherever necessary, but 
the physical background for the mathematical expressions is always emphasized. From 
the mathematician’s point of view such treatment is not elegant, but it is much to this 
reviewer’s liking. 

This book is particularly recommended to all technical men recently employed in 
geophysical work. It should also prove useful to geologists associated with geophysical 
surveys. 

To summarize: This is a good book, well written, and well illustrated. 

W. T. Born 


SUMMARY 


It is noteworthy fact that three very excellent books dealing with geophysical pros- 
pecting should appear in the United States almost simultaneously. A great deal of 
literature in the field of applied geophysics has been accumulating scattered through 
the physical, geological, mining and geophysical journals and the books which have been 
available up to now really belong to an earlier era of this rapidly advancing art. All 
three authors have substantial practical experience in geophysical prospecting. All of 
them have taught the subject to college students and all of them have made contribu- 
tions to its technical literature. Nettleton’s book, reviewed here by Dr. Born, restricts 
itself to prospecting for oil and is, if that restriction is acceptable, the best introductory 
text on the subject now available. At the same time the book is at least equally interest- 
ing to the specialist, especially his treatment of the gravity method of prospecting, of 
which he is an outstanding practitioner. 

As is perhaps to be expected, all the books strongly reflect the known special in- 
terests of their respective authors. This is perhaps best illustrated by the following 
table which illustrates the proportions of space allocated to major subjects: 











Magnetic Gravity Seismic Electrical* 
Jakosky 12.2% 12.7% 24.9% 24.7% 
Heiland 14.2% 23.3% 18.1% 20.4% 
Nettleton 16.9% 32.7% 30.0% 5.0% 


* Exclusive of well logging. 


The seismic method is at present, has been and at least so long as oil is the pre- 
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dominant objective of geophysical prospecting, is likely to be the most widely practiced 
of geophysical methods. If degree of utility is the criterion for space allocation, the 
seismic method undoubtedly deserves first place. As a matter of history, the gravity 
method unquestionably ranks second on the same basis. Nettleton’s 5% of space to 
electrical prospecting (exclusive of well logging) is more in line with the limited aggre- 
gate of such work done in prospecting for oil. Jakosky’s and Heiland’s substantially 
greater space allocation recognizes the great volume of published literature and of 
patents on electrical prospecting. 

It seems unfortunate that the titles of Jakosky’s and Heiland’s books are not inter- 
changed because, although the difference is not very great, Jakosky’s subject seems to 
the reviewer to be geophysical exploration and Heiland’s, exploration geophysics. In 
each book a geophysicist may expect to find things he may want and which he will not 
find in the others and probably he can not afford to leave any one of them off his book- 


shelf. 
E. A. ECKHARDT 


Lecture Notes on Practical Petroleum Geophysics, by Garrett Kemp. 


This small booklet contains lecture notes from a ‘“‘Conference on Practical Petro- 
leum Geophysics,” held at Texas A & M College, April 4-6, 1940. 

One gathers the impression from reading these notes that the Conference was quite 
interesting. There are some forty-five well chosen and informative figures. Unfortunate- 
ly, however, the authorities, who were responsible for the transfer of this material from 
the spoken lectures to the printed booklet, failed. As a result, so many errors have crept 
in that no student should read it until he is very well informed about practical geophys- 
ics. 

The recent publication of three excellent books on practical geophysics makes it 
unnecessary to endeavor to point out and correct each error in this booklet. Moreover 
this would require a review practically as long as the original. Therefore only one typical 
example will be cited from each section. 

Considerable space is devoted to a comparison of the torsion balance, gravimeter 
and pendulum. First, relative sensitivities are compared. The torsion balance measures 
the gravitational gradient, the unit of which (page 21) ‘‘is named for Eétvés and is 
0.000,000,001 dyne per centimeter. It will be seen that the Eétvés unit is 100,000 times 
smaller than the milligal which is the unit used in pendulum and gravimeter work.” 
Consequently (page 22), ‘‘it would be unscientific to compare the practical values of 
these instruments with the torsion balance . . . simply because the pendulum and the 
gravimeter must have their respective sensibilities increased at least ten thousand times 
which would be practically impossible.” As a result (page 27), ‘‘structural features can 
be located with the torsion balance that could not be even suggested by the present 
gravimeter because the sensitivity is deficient from one hundred to one thousand times 
that of the torsion balance.’’ Because they failed to comprehend this fact ‘‘certain 
pseudo-geophysicists were misapplying torsion balance gradients in such a way as to 
create the impression that after all torsion balance surveys are identical with those of 
of the gravimeter” (page 18), that is by constructing isogams. The seriousness of their 
error is made clear bya numerical example (page 25) where it is shown that a gravim- 
eter would give 0.012 as the difference between the gravitational forces at two hy- 
pothetical points whereas from torsion balance data it can be seen that the increase in 
gravitational force between these two points is 0.000,549. 
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The pendulum is more valuable than the gravimeter because (pages 12, 19, 16) it 
gives the absolute value of gravity. Since (page 16) “‘the cost of making a survey is 
dependent upon the results and the financial returns,” is not surprising that (page 15) 
“the pendulum equipment is much more costly than any other gravitational method 
because the investment in equipment is greater.” 

The magnetic section contains fewer errors than the rest, but even in it there are 
such statements as (page 36) “‘the magnetic susceptibility of a vacuum is taken as 
unity. All organic compounds measured have values less than unity, that is negative.” 

The section on temperature includes a rather confusing paragraph (page 48). ‘This 
work suggests that oil well temperatures conform to a simple law, namely, the tempera- 
ture at the depth, D, is equal to the computed mean annual temperature for the area 
plus the product of the increase in temperature per foot and the depth in feet. However, 
the surface temperature can be used for all practical purposes since the temperatures 
used in computation can be taken at small intervals from the surface, some of which 
will indicate a constant temperature, independent of surface temperature changes 
accurate enough in the study of well temperatures in a particular area.” 

Most seismic prospectors would doubt that (page 52) ‘‘Seismograph work would be 
simplified and made more dependable if an explosion could be devised so as to produce 
a series, instead of an impulse, at equal intervals of time then the average velocity 
could be measured on the film record.”’ Nor is it usually believed (pages 53 and 55) that 
first arrivals are always waves which have travelled along shallow horizontal layers. 

In the section on electrical methods, methods depending on natural currents (page 
60) are said to be valueless and yet those depending on natural potentials (page 62) 
have ‘‘very favorable possibilities.” The Wenner-Gish method (page 63) “‘is limited as 
to the depth of penetration and sufficient sensitivity to differentiate between geologic 
strata,” and consequently is of little practical value, but (page 65) the modification 
known as Lee’s Partitioning Method “‘is expected to take its place as the most promis- 
ing geophysical method that has been developed without exception to any other geo- 
physical method. It is stated that (page 65) ‘“The following is a list of accomplished 
facts, possible with this method, based upon actual field work. 

(1) Geologic strata or beds mapped accurately by determining characteristic 
changes in resistivity at depths measured accurately. 

(2) Geologic structure accurately mapped. In some places to 7000 ft. depth. 

(3) Changes in impregnation (or related porosity) or strata mapped accurately. 

(4) Faults in any kind of strata mapped directly and accurately. 

(5) Gas located, depth and areal limits mapped. 

(6) Salt water level determined accurately. 

(7) Salt water contacts with oil, gas, fresh water, sand, shales, rocks, & etc. de- 
termined and mapped accurately. 

The Lee Partitioning Method is capable of obtaining more geological information 
than is possible by any other geophysical method at no greater cost.” It would be 
interesting to know more about this ‘‘actual field work.” 

W. M. Rust, Jr. 


Strength and Structure of the Earth, by Reginald Aldworth Daly. Prentice-Hall Geology 
Series, New York 1940. 434-+1X pp. $3.50. 


For the most part, the contents of this book are not what one would expect from 
its title. The only part that bears directly on the subject of the title is the last one of 
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the twelve chapters. The remainder is devoted to the assembly of the information on 
which the estimation of the strength of the earth’s outer crust depends. Since such 
estimates are based almost entirely on the degree of departure from isostasy, the book 
is very largely devoted to the measurements of gravity and of deflections of the vertical 
and their implications in terms of the shape of the earth and the degree to which topog- 
raphy is balanced by density variations within the crust. 

It appears that Professor Daly has made a very thorough study of the original 
sources of the fundamental data on which the determination of the earth’s strength 
depend and has then assembled the material into this volume. The scope of the material 
thus assembled is indicated by the list of special acknowledgments to institutions and 
authors: this includes the United States Coast & Geodetic Survey, the Canadian Do- 
minion Observatory, the Survey of India, the Finnish Geodetic Institute, and the Baltic 
Geodetic Commission, and Messrs. Bowie, Cassinis, Glennie, Gutenberg, Heiskanen, 
Jeffreys, Jung, Lambert, Matuyama, Sauramo, Vening Meinesz, and others. 

The book includes many details of gravity surveys throughout the world and their 
discussion in tems of different schemes of isostasy. The presentation includes a great 
many quotations from the acknowledgments listed above and they constitute some 
20% of the total content of the book. The wide range of sources covered serves to com- 
pare the details of the various systems of isostasy which have come about through in- 
vestigations by different men of data from widely different parts of the earth’s surface. 
In particular, those who are familiar with the complete zone-by-zone isostasy of Hay- 
ford and Bowie will come to realize that this is probably far from a physical fact and 
is largely a fiction of mathematical convenience. A much broader system of equilibrium, 
which is combination of the Pratt and Airy systems of isostasy, apparently is nearer to 
reality. That this feeling is developed after reading the book is probably a reflection of 
many quotations from Heiskanen and Professor Daly’s evident belief that Heiskanen’s 
ideas, including the tri-axial ellipsoid, are the nearest current approach to a true rep- 
resentation of conditions in the crust. 

Considerable attention is given to the variation in the isostatic anomalies depend- 
ing upon the variations of the systems of isostasy on which the calculations are based 
and upon the form of the geoid to which the anomalies are referred. An attempt is made 
to reconcile some of the apparent inconsistencies in the strength requirements for the 
earth’s shells by consideration of the different bases on which these anomalies have 
been computed in the past and reference of all to a more uniform system. Thus it is 
shown that the apparent weakness of the earth’s crust and high fluidity of the under- 
lying heavy asthenosphere, as indicated by the recovery of Fennoscandia and Labrador 
from the burden of the Pleistocene ice cap in contrast with the apparent strength of the 
crust, as indicated by the large anomalies in peninsular India, are not grossly inconsis- 
tent when both sets of anomalies are computed from the Heiskanen tri-axial ellipsoid 
and when consideration is given to the different areas of the anomalies in the two cases 
and the consequent demands upon the strength of the crust supporting them. 

As a review of the gravitational tests of isostasy and collection of data from all over 
the world, the present book may be considered somewhat as a sequel to Bowie’s Jsos- 
tasy. However, Professor Daly’s work is broader in that it gives rather careful and 
quite unbiased consideration to different investigators with conflicting points of view 
and conclusions. The almost complete absence of mathematics may appeal to the geo- 
logical or general reader, but in some cases attempts to explain mathematical equations 
in words are either inaccurate or inadequate to the extent that they might better have 
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been omitted. This is particularly true with respect to the formula for gravity on the 
surface on the earth (p. 109) and in the use of the Clausius—Clapeyron equation (p. 
412). There are a few minor quantitative inaccuracies and errors of description, partic- 
ularly of pendulum apparatus (p. 112), and some places in which quotations from other 
writers give an incorrect impression when set out by themselves. However, Professor 
Daly has done a real service in bringing together a critical review of the present status 
of isotasy and its bearing on the strength of the earth’s crust and in showing the limits 
to the depths and properties of the outer part of the earth which are derivable from 


present gravity and deflection observations. 
L. L. NETTLETON 


Annual Reviews of Petroleum Technology, Vol. 5 (covering 1939), published by the Insti- 
tute of Petroleum, c/o the University, Edgbaston, Birmingham 15, England. 444 
pp., U.S.A. $2.70 incl. postage. 


It is of interest that the present volume, the fifth in the series of annual reviews 
of developments in the science and technology of petroleum, has been published under 
the stress of war conditions in England. Inevitable restrictions and voluntary censor- 
ship has led to omissions by the authors of the twenty-eight chapters comprising the 
volume. However, the reader will find a comprehensive review covering all phases of 
the industry from, regional exploration and development, drilling, production, 
‘through refining and including engine developments. During 1939 progress was prob- 
ably most marked on the refining side of the industry where new processes of polymer- 
ization, alkylation, cracking and reforming rapidly developed. Thus a considerable 
portion of the review covers such developments. 

Geophysics is covered in a brief six pages by J. H. Jones. Hence the volume is not 
recommended to geophysicists who are interested only in their own field. On the other 
hand to one interested in following the development of other phases of the industry, 
it provides a concise résumé of progress and modern trends presented by some forty 
competent reviewers. Also of interest to students of the industry are the lists of refer- 
ences following each chapter comprising not less than 2700 items. 

R. D. WycKoFF 








PATENTS* 
ELECTRICAL PROSPECTING 


U.S. No. 2,217,780, T. S. West and C. C. Beacham, Iss. 10/15/40, App.. 7/26/37. . 


Method of Electrical Geophysical Prospecting and Apparatus for Practicing Said 
Method—Distribution of current through the earth is changed while recordings are 
made of potential differences picked up by several potential electrodes. Vacuum tube 
potential indicators can be set to zero deflection so that changes in potentials can be 
observed directly. 


GRAVIMETRIC PROSPECTING 


U.S. No. 2,217,123, J. D. Malmqvist, Iss. 10/8/40, Appl. 10/11/38, Assign. Bolidens 
Gruvaktiebolag, Sweden. 


Apparatus for Determing the Alterations of the Horizontal Components of the Force 
of Gravity—The angle between the suspension wires of two spaced plummets is deter- 
mined by measuring the electrostatic capacity between the plummets and adjoining 
fixed plates. 


U.S. No. 2,217,361, H. M. Evjen and D. S. Muzzey Jr., Iss. 10/8/40 Appl. 5/21/38, 
Assign. Shell Development Company. 


Gravity Measuring Instrument—A gas-spring gravimeter of the general astatized 
type in which the system is balanced on a knife-edge allowing it to tilt as the liquid is 
displaced by changes of gravity. Tilting increases the normal displacement of liquid 
thus increasing sensitivity. Stability may be regulated by shifting a weight to raise 
or lower the center of gravity. 


U.S. No. 2,218,140, F. E. Wright and J. L. England, Iss. 10/15/40, Appl. 12/15/38, 

Assign. Carnegie Institution of Washington. 

Gravity Torsion Balance—A weighted arm radially supported at the center of a 
horizontal helix is deflected to a definite position by rotating the helix-supporting frame 
a measurable amount about the helix axis to determine the force of gravity. Tempera- 
ture control, compensation and indication are provided. Clamping means and level 
viewing means are included. 


U.S. No. 2,220,199, J. L. Bible and R. H. Ray, Iss. 11/5/40, Appl. 12/27/38, Assign. 
Stanolind Oil and Gas Company. 


Clamping Device for Force Responsive Element—Bar-type bifilar gravimeter is 
clamped between two pairs of faces which approach from opposite sides to clamp the 
bar in approximate operating position. 


U.S. No. 2,221,480, G. A. Ising, Iss. 11/12/40, Appl. 2/11/38. 


Gravity Measurements—The support of an astatized pendulum is periodically tilted 
a fixed amount while deflections are recorded. The degree of astatization is indicated 
by the amount of shift of the pendulum. 


*Quarterly listings provided by Gary Muffly, G. R. & D. Co. 
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SEISMOGRAPH PROSPECTING 
U.S. No. 2,215,297, J. E. Owen, Iss. 9/17/40, Appl. 2/16/40, Assign. Geophysical Re- 
search Corporation. 
Seismic Surveying—Surface waves are cancelled out by combining outputs of two 
geophones placed at different angles with respect to vertical. Phase and magnitude of 
of the surface waves received are adjusted for cancellation. 


U.S. No. 2,216,452, J. E. Owen, Iss. 10/1/40, Appl. 5/27/38, Assign. Geophysical Re- 
search Corporation. 
Seismic Surveying—Method of cancelling surface waves by phase and amplitude 
adjustments of outputs of horizontal and vertical seismometers. 


U.S. No. 2,217,720, J. F. Anderson, Iss. 10/15/40, Appl. 9/19/38, Assign. Sun Oil 
Company. 
Apparatus for Solving Seismographic Problems—Reflection problems are solved 
with a mechanical device which locates the reflecting surfaces geometrically. 


U.S. No. 2,217,806, Gary Muffly, Iss. 10/15/40, Appl. 7/29/37, Assign. Gulf Research 
& Development Company. 
Seismograph Pros pecting—Tuning of a seismograph amplifier is varied as a func- 
tion of time. 


U.S. No. 2,217,828, R. D. Wyckoff, Iss. 10/15/40, Appl. 7/16/37, Assign. Gulf Re- 
search & Development Company. 


Seismograph Prospecting A pparatus—A heterodyne type amplifier converts seismic 
signals to high frequency side bands. A sharply tuned filter selects the frequencies in 
one side band which correspond to the desired seismic frequencies and a detector con- 
verts them back to the original low frequency. The filter frequency is fixed and identical 
in all channels and tuning is accomplished in all channels simultaneously by changing 
the frequency of the heterodyne oscillator common to all. 


U.S. No. 2,219,508, L. F. Athy and E. V. McCollum, Iss. 10/29/40, Appl. 8/31/38, 

Assign. Continental Oil Company. 

Method and Apparatus for Seismic Surveying—Seismograph reflection setups have 
one line or two intersecting lines of detectors offset from the shot point by a distance of 
75 to 1000 feet, thus providing continuous control on the reflecting horizon but remov- 
ing all detectors from the immediate vicinity of the shot and the disturbing ground- 
roll zone. 

GEOCHEMICAL PROSPECTING 
U.S. No. 2,223,183, FE. B. Peck, Iss. 11/26/40, Appl. 9/2/37, Assign. Standard Oil 

Development Company. 

Prospecting for Oil—Soil gas samples are burned and the temperature rises due to 
their combustion are measured. Samples showing greatest combustion effect are tested 
for carbon-hydrogen ratio by correlating heat of combustion with CO; content. 


U.S. No. 2,223,785, G. L. Hassler, Iss. 12/3/40, Appl. 12/20/38, Assign. Shell Develop- 
ment Co. 
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Gas Sampling System—A container for storing gaseous mixtures designed for de- 
livering such mixtures to and from the containers without contamination or dilution 
by atmospheric air or impurities. Comprises two concentric gas-tight envelopes at least 
the inner one being made of flexible material whereby suitable manipulation of gaseous 
pressure permits filling or emptying either container. 


WELL LOGGING 
U.S. No. 2,212,273, Oscar Martienssen, Iss. 8/20/40, Appl. 10/12/38. 
Arrangement for Measuring the Local Specific Resistance of Borehole Strata—Well 
logging electrode consists of a heavy lead weight which will make contact with the bore 


wall. Weight is covered with insulating jacket except for a small central section so that 
contact is localized. 


U.S. No. 2,212,274, Oscar Martienssen, Iss. 8/20/40, Appl. 10/21/38. 


Method of Exploring the Porosity of Geologic Strata Traversed by Boreholes—Self 
emf’s due to porosity are logged by measuring the current flow while a small D.C. po- 
tential is applied. Counter emf. from induced electro-osmotic flow causes current to 
decrease opposite porous formations. 


U.S. No. 2,214,674, J. T. Hayward, Iss. 9/10/40, Appl. 1/29/38. 


Method of Logging Wells—A well is logged during drilling by analysis of the drill 
mud before and after it passes through the well. Rate of mud flow, depth, etc., are ob- 
served so that a portion of the fluid entering the well can be identified after returning 
to the surface, and depth at which it picked up disintegrated rock may be determined. 
Viscosity, gas, oil or salt content, etc., of samples may be studied. Mud pressure is kept 
higher than formation pressure. 


U.S. No. 2,214,786, B. P. Bishop, Iss. 9/17/40, Appl. 8/15/38. 


Apparatus for Logging Holes while Drilling—Well is logged while being drilled by 
measuring self e.m.f. between uninsulated drill stem and a surface electrode. 


U.S. No. 2,216,340, R. D. Elliott, Iss. 10/1/40, Appl. 6/25/37, Assign. Lane-Wells 
Company. 
Electric Cable—Well logging cable. Each conductor has a stranded steel sheath for 
magnetic shielding. 


U.S. No. 2,217,768, J. M. Pearson, Iss. 10/15/40, Appl. 11/5/37, Assign. Sperry-Sun 

Well Surveying Company. 

Apparatus for Orienting Cores—A core from a well bore is tested by rotating it in 
the vicinity of an astatic magnetometer while recording deflections. The core is then 
reversed to obtain a second record which is compared to the first record to separate 
polarization and susceptibility effects. 


U.S. No. 2,219,273, S. A. Scherbatskoy, Iss. 10/22/40, Appl. 6/16/39, Assign. Well 
Surveys, Inc. 
Well Logging by Measurements of Radioactivity—Radioactivity well logging meas- 
urements are made with an ionization chamber whose output is commutated, ampli- 
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fied with an A.C. amplifier and recorded at the surface. Ionization chamber uses 
nitrogen under about 20 atmospheres of pressure. 


U.S. No. 2,219,274, S. A. Scherbatskoy, Iss. 10/22/40, Appl. 10/16/39, Assign. Well 

Surveys, Inc. 

Well Survey Method and A pparatus—A null method of recording wherein the D.C. 
output of an ionization chamber in a well is balanced out by remote control from the 
surface and any unbalance is indicated by means of an A.C. amplifier with commutated 
D.C. input. 


U.S. No. 2,220,070, C. B. Aiken, Iss. 11/5/50, Appl. 3/24/36, Assign. Schlumberger 
Well Surveying Corp. 
Method of and Apparatus for Magnetically Exploring Earth Strata—Inductive 

coupling between two or more coils is indicated while the coils are moved through a 

well in spaced relationship to log the well. Setup may be balanced or null-reading, and 

may indicate phase as well as amplitude. Frequency variations of an oscillator in the 
well may be measured at the surface. 


U.S. No. 2,220,205, S. E. Buckley, Iss. 11/5/40, Appl. 3/31/39, Assign. Standard Oil 
Development Company. 
Method of Locating Detectable Cement in a Borehole—Cement ina bore hole is located 
by adding radioactive material to the cement and logging the well for emitted rays. 


U.S. No. 2,220,509, Folkert Brons, Iss. 11/5/40, Appl. 1/30/40, Assign. Shell Develop- 
ment Company. 
Process and Apparatus for Exploring Geological Strata—Well is logged by passing 
a neutron source through it while detecting neutrons diffused by the surrounding forma- 
tions. Returning neutrons disintegrate a selected element such as cadmium, and the 
disintegration products are detected with a Geiger counter or ionization chamber. 


U.S. No. 2,220,788, R. W. Lohman, Iss. 11/5/40, Appl. 6/24/35, Assign. Schlumberger 

Well Surveying Corp. 

Method and Apparatus for Investigating Subterranean Strata by Means of Electro- 
magnetic Measurements—Well is logged for magnetic or inductive effects. Induction 
between exploring coils traversing the bore hole is studied with respect to amplitude, 
phase, or distortion. Earth’s field may also be logged. 


U.S. No. 2,221,951, W. D. Mounce and W. M. Rust, Jr., Iss. 11/19/40, Appl. 12/24/37, 
Assign. Standard Oil Development Company. 


Making Electrical Measwrements—A transformer placed near logging electrodes in 
a well matches the low electrode impedance to the high measuring-circuit impedance 
so that changing cable inductance will not affect measurements. A condenser is pro- 
vided to allow simultaneous measurement of natural D.C. earth potentials. 


U.S. No. 2,222,136, C. B. Bazzoni and J. Razek, Iss. 11/19/40, Appl. 8/2/38, Assign. 
Sperry-Sun Well Surveying Company. 


Electrical Prospecting Method and A pparatus—A well logger records variations as a 
variable-density trace which is reproduced later in the form of a curve. Variations 
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may be balanced against a simultaneously recorded standard-exposure trace to cancel 
spurious effects of variable light source or film speed. 


U.S. No. 2,222,149, E. Lipson, Iss. 11/19/40, Appl. 11/28/39, Assign. Casolog, Inc. 


Method and A pparatus for Logging Wells—Electrical well logging setup uses a single 
well electrode and two surface electrodes connected in a bridge circuit. 


U.S. No. 2,222,182, W. D. Mounce and W. M. Rust, Jr., Iss. 11/19/40, Appl. 8/21/37, 
Assign. Standard Oil Development Company. 
Earth Impedance Measuring Device—Well logger or the like has two electrodes con- 
nected in the feedback circuit or in a tuned circuit of a vacuum tube oscillator. Changes 
of impedance vary the strength of oscillation which is indicated by a plate current meter 


at the surface. 


U.S. No. 2,222,608, R. D. Elliott, Iss. 11/26/40, Appl. 5/17/37, Assign. Schlumberger 

Well Surveying Corp. 

Electrical Surveying in Drill Holes—Electrical well-logging circuit with compensa- 
tion for drift of the recording meter due to changes in electrical characteristics of the 
cable as it is unwound. Current and potential circuits may each use one long electrode 
and one short electrode. Cable shield may be used as one electrode. Shields may be used 
between conductors. Electrode assembly details are described and claimed. 
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ANNOUNCEMENT OF ANNUAL MEETING—1941 


The Annual Meeting of the Society of Exploration Geophysicists will be held at 
the Rice Hotel, Houston, Texas, April 1 to 3, inclusive. 


TENTATIVE PROGRAM 


Tuesday, April 1—A.M. Business Meeting and Technical Session. 
P.M. Technical Session. 
Wednesday, April 2—A.M. Technical Session. 
P.M. Technical Session. 

Thursday, April 3—A.M. Joint Session with A.A.P.G. and S.E.P.M. at which time the 
presidential addresses and other papers of general interest 
will be presented. 

P.M. No Technical Session in order that members may hear pa- 
pers of general interest on the A.A.P.G. program. 


PROGRAM COMMITTEE 


E. A. Eckhardt— _ Atlantic Coast Area 
J. A. Sharpe— Tulsa Area 
C. A. Heiland— Rocky Mountain Area 


Henry C. Cortes— _  Dallas-Ft. Worth Area 
Herbert Hoover, Jr.—Pacific Coast Area 
Frank Goldstone— Houston Area 


T. I. Harkins— Houston Area 
W.M. Rust, Jr— Houston Area 
H. B. Peacock— Houston Area (General Chairman) 


IT IS URGED THAT ALL MEMBERS DESIRING TO PRESENT PAPERS 


NOTIFY ANY MEMBER OF THE PROGRAM COMMITTEE AS SOON AS 
POSSIBLE. 
HOTEL ACCOMMODATIONS 
About 1500 rooms have been reserved at the various hotels of the city for the com- 
bined conventions. However, it is necessary that individual reservations be made before 
March 20. After that date all unreserved rooms will be made available to the general 
public. The following hotel information is given for your convenience: 


Minimum Rates 


No. of No. Blocks Without Bath With Bath 

Hotel Rooms _ from Rice Single Double Single Double 
Rice 1000 $2.50 $4.00 
Lamar 500 . 2.50 4.00 
Texas State 400 3 2.50 4.00 
San Jacinto 300 3 $1.50 $2.50 2.50 4.00 
Ben Milam 250 6 2.50 4.00 
Sam Houston 200 2 2.50 4.00 
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Minimum Rates 


No. of No. Blocks Without Bath With Bath 

Hotel Rooms from Rice Single Double Single Double 
William Penn 175 C 2.50 4.00 
Cotton 175 3 2.50 4.00 
Auditorium 175 2 2.50 4.00 
Warwick 153 50 2.50 4.00 
De George 150 7 $2.00 $3.00 2.50 4.00 
Plaza 102 50 2.50 4.00 


MEMBERSHIP APPLICATIONS APPROVED FOR PUBLICATION 


The Executive Committee has approved for publication the names of the following 
candidates for membership in the Society. This publication does not constitute an 
election, but places the names before the membership at large. If any member has 
information bearing on the qualifications of these nominees, he should send it to the 
Secretary within thirty days. (Names of sponsors are placed beneath the name of each 
nominee.) 


ACTIVE 


Dorsey Augustus DeWoody 

E. A. Eckhardt, R. D. Wyckoff, L. W. Gardner 
Jardel Borges Ferreira 

Herbert Hoover, Jr., Wayne H. Denning (Article III-C-1) 
Holloway Halstead Frost 

L. M. Mott-Smith, E. V. McCollum, G. P. Wilson 
A B Hamil 

O. S. Petty, R. B. Pratt, W. Harlan Taylor 
John Brackett Hersey 

Maurice Ewing, W. B. Agocs (Article III-C-1) 
Paul Hugh Jeffers 

Arland I. Innes, Andrew Gilmour, L. Y. Faust 
Lewis Owen Kelsey 

W. R. Ransone, Earl Thomas, C. V. A. Pittman 
Noyes Darling Smith 

J. P. Woods, H. M. Evjen, D. S. Muzzey, Jr. 
Richard Hanawalt Tallman 

A. B. Bryan, George E. Wagoner, Robert L. Craig 
Garvin Lawrence Taylor 

Daniel Silverman, Joseph A. Sharpe, R. Clare Coffin 
Eugene England Wilson 

James L. Morris, W. H. Schley, W. W. Newton 


ASSOCIATE 


Charles Carpenter Bates 

George E. Wagoner, A. B. Bryan, R. W. Gemmer 
Horace Chesterfield Davis 

G. R. Brotherhood, Ray K. Carter (Article ITI-C-1) 
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Burleson Graham 

R. A. Pratt, R. A. Weisbrich, T. L. Allen 
C. Russell Hale 

Wayne H. Denning, A. M. Eichelberger, Jr., Herbert Hoover, Jr. 
Gerald Willis Hostetler 

T. L. Allen, John F. Imle, W. O. Bazhaw 
John Byron Jones 

G. F. Kaufmann, C. A. Heiland (Article ITI-C-1) 
Lynn Alton Smitherman 

W. Lee Moore, W. O. Bazhaw, T. L. Allen 
Charles Everett Thompson 

W. O. Bazhaw, R. B. Pratt, R. A. Weisbrich 
Wilber H. Young, Jr. 

Sylvain Pirson, W. R. Ransone, Ray K. Carter 


TRANSFER TO ACTIVE 
(Executive Order G-11) 
Willis Oliver Bazhaw 
James Dungan Berwick 
Raiford H. Burton 
James Edward Gunn 
E. Darrell Williams 
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The ballot reprinted below was mailed all active members of the Society of Explora- 
tion Geophysicists on December 15, 1940. Active members who failed to receive bal- 
lots may secure them by notifying J. F. Gallie, Business Manager, P. O. Box 777, 
Austin, Texas. Only active members who have paid 1941 dues are eligible to vote. 


OFFICIAL BALLOT 


FOR THE YEAR MARCH, 1941—MARCH, 1942 
For the Election of Officers for the 
SOCIETY OF EXPLORATION GEOPHYSICISTS 








FOR PRESIDENT 
Vote Once Only 





























H. B. PEACOCK O 
ROLAND F. BEERS C) 
O 
FOR VICE-PRESIDENT 

Vote Once Only 
FRANK GOLDSTONE O) 
HENRY SALVATORI oO 
O 

FOR SECRETARY-TREASURER 

Vote Once Only 
WILLIAM M. RUST, JR. O 
JOSEPH A. SHARPE O 
O 

FOR EDITOR 

Vote Once Only 
RICHARD L. DOAN 
R. D. WYCKOFF CJ 
L 











After voting, place in accompanying official envelope, which should then be signed and 
mailed to the 


SECRETARY-TREASURER, 
ANDREW GILMOUR 
P. O. Box 2040 
Tulsa, Oklahoma 



























PERSONAL ITEMS 


H. W. Srratey III, has the distinction of having placed all of his senior class 
(June, 1040) in permanent positions within the profession by July 1 and his entire junior 
(June, 1941) class in temporary summer employment. The enrollment in the depart- 
ment of geology at Baylor has increased 39% since Mr. Straley assumed its direction, 
and the number of B.S. majors has doubled. 


Advice has been received that WARREN A. ALEXANDER is now in Ploesti, Roumania, 
with Romano-Americana. 


Howarp L. Coss, newly-elected member, may be addressed at Box 491, Hanover, 
New Hampshire. He is in attendance at Dartmouth College. 


GEORG MARKSTEIN is no longer connected with Compania de Petroleo Shell de 
Colombia. He should be addressed at Calle 35, No. 17/21, Bogota, Colombia. 


GirForD E. WHITE has advised the Society of his new address—273 Beacon 
Street, Boston, Massachusetts. 


HENRY SCHOELLHORN III, a recent graduate of the Colorado School of Mines, 
has accepted employment with the Gulf Research and Development Company. He 
may be addressed c/o Seismic Party #20, P. O. Drawer 552, Liberty, Texas. 


C. T. MAcALLISTER is now with the Independent Exploration Company, 2011 
Esperson Building, Houston, Texas. 


Matvin G. HorrMan, whose ‘‘Structural History of Billings Field, Noble County, 
Oklahoma, Interpreted in Terms of Isostasy” appeared in the November A.A.P.G. 
Bulletin, is geologist with Geochemical Service Corporation, 321 S. Detroit, Tulsa, 
Oklahoma. 


The Petroleum Engineer has announced the election of the following officers: W.L. 
Love, President; K. C. ScLaTEr, First Vice-President and Editor; T. J. CROWLEy, 
Vice-President and Advertising Manager; and W. T. Bryan, General Manager and 
Secretary-Treasurer. 


Roy LEE GALLaway, recently returned from work in Egypt, and employed by 
the Texas Company, is temporarily located in Alexandria, Louisiana. 


T. C. RicHarps, in Burma for some time past, has returned to England, where he 
may be addressed at 104, Archery Road, Eltham, London. 


Francis H. Capy, of the Carter Oil Company, is in Aberdeen, Mississippi, with 
his postal address being Box 252. 


GARDNER P. WILSON, Junior Detailer with the Lockheed Aircraft Corporation of 
Burbank, California, is located at 946 N. Catalina in that city. 


VINCENT EMANUEL, in Italy for the past few months, is returning shortly to the 
United States. 


R. G. PATERSON, recently-elected member, is a Lieutenant with the R. C. N. V. R. 
and may be reached through the Fleet Mail Office, Halifax, Nova Scotia, Canada. 
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Wit1aM B. Acocs has left the employ of the Gulf Research and Development 
Company to do graduate study in Physics at Lehigh University. His address is 940 
Broadway, Bethlehem, Pennsylvania. 


The Jackson, Mississippi, office of Geophysical Service, Inc., has been moved to 
207 100 East Pearl Building, Ewrn D. Gasy reports. 


FRANK F. LEHMAN is working in Barcelona, Venezuela, for the Mene Grande Oil 
Company, Apartado 45. 


F. H. AGEE may be addressed at P. O. Box 413, Longview, Texas. 


H. RUTHERFORD is now located in Centralia, Illinois, where his mailing address is 
P. O. Box 96. 


Ray K. Carter has returned from South America to his position with Geophysical 
Service, Inc., 1311 Republic Bank Building, Dallas, Texas. 


E. W. FISHER is serving as an Electrical Engineer with the War Department. His 
local address is 200 Courtland Avenue, Topeka, Kansas. 


Pau H. Boots is employed by the Mene Grande Oil Company, Seismograph 
Division, Apartado 45, Barcelona, Venezuela. 


E. R. SELBy has accepted employment as an Engineer with the Magnetest Corpora- 
tion, 3504 Atlantic Avenue, Long Beach, California. 


Epcar J. STULKEN is acting as an Instructor in Mathematics at Allen Academy, 
Bryan, Texas. 


T. De PiyPER has left Burma for Pangkalan Brandan, Sumatra, where he is work- 
ing for N. V. Bataafsche Petroleum Maatschappij. 


W. B. WILEy, of the Gulf Research and Development Company, has been trans- 
ferred from Watonga, Oklahoma, to Columbia, Mississippi. 


A. A. Hotarp is located in Houma, Louisiana, at 352 Ruth Street. 


A. N. Esy, of the Superior Oil and Gas Company, is temporarily located in Colum- 
bus, Texas. 


J. W. Tuomas, Jr., of the Stanolind Oil and Gas Company, has been transferred 
from Canadian, Texas, to Garden City, Kansas. 


L. W. Brau, J. P. ScoumacHEeR, PAUL WEAVER, FREDERICK ROMBERG, W. O. 
BazHaw, J. BRIAN Esy, F. W. RotsHausEN, E. De Gotyer, R. B. Cort, F. Juttus 
Fous, Bruce H. Hartton, Cuas. Girt Moraan, H. W. Stratey III, Jostan TAyYLor, 
LEO W. Kowz, C. V. A. Pitrman, and Marion J. Moore were among Society members 
attending the 53rd Annual Meeting of the Geological Society of America, which was 
held December 26-28 in Austin, Texas. 


RoseErt G. NISLE is now associated with the Dictograph Sales Corporation, 408 
West Building, Houston, Texas, as Telephone Engineer. 


R. B. Cort, of Geophysical Service, Inc., has returned to the United States from 
the Netherlands East Indies. 


Dwicut E. Warp, of the Carter Oil Company, has been transferred from Phila- 
delphia, Mississippi, to Phillipsburg, Kansas, where his postal address is Box C. 
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GrEorGcE B. KalsER is now serving as a Captain at Fort Monmouth, Red Bank, 
New Jersey. 


Homer Roserts, of the Petty Geophysical Engineering Company, has returned 
from Burma to San Antonio, Texas. 


R. A. Watson, of the National Geophysical Company, has been transferred from 
Alexandria to Winnfield, Louisiana, with the postal address Box 407. 


J. C. Horrman, of the United Geophysical Company, is located at 317 Haberfelde 
Building, Bakersfield, California. 


NErL W. Manv, of Geophysical Service, Inc., has returned from Bahrein Island to 
Dallas, Texas. 


SHERWIN F. KELLy, who recently gave several lectures and field demonstration 
in exploration geophysics to students of the Montana School of Mines in Butte, haS 
been named Walker-Ames Lecturer in Geophysics at the University of Washington 
in Seattle. During his visiting lectureship he will give eight lectures on the subject, 
emphasizing the practical importance of the geophysical sciences and the interrelation- 
ship with geology. 


see OF THE OWNERSHIP, MANAGEMENT, CIRCULATION, ETC., 
REQUIRED BY THE ACTS OF CONGRESS OF AUGUST 24, 1912, 
AND MARCH 3, 1933 


of GEOPHYSICS, published quarterly at Menasha, Wisconsin, for October 1940. 


State of Texas 
County of Travis 

Before me, a Notary Public in and for the State and county aforesaid, personally appeared 
J. F. Gallie, who, having been duly sworn according to law, deposes and says that he is the Business 
Manager of the Geophysics and that the following is, to the best of his knowledge and belief, 
a true statement of the ownership, management (and if a daily paper, the circulation), etc., of the 
aforesaid publication for the date shown in the above caption, required by the Act of August 24, 
1912, as amended by the Act of March 3, 1933, embodied in section 537. Postal Laws and Regula- 
or printed on the reverse of this form, to wit: 

. That the names and addresses of the publisher, editor, managing editor, and business 
ianiaaea are Publisher, George Banta Publishing Company, Menasha, Wisconsin; Editor, 
Wyckoff, Box 2038, Pittsburgh, Pennsylvania; Managing Editor, R. D. Wyckoff, Box 2038, 
Pittsburgh, Pennsylvania; Business Manager, J. F. Gallie, Box 777, Austin, Texas. 

2. That the owner is (If owned by a corporation, its name and address must be stated and 
also immediately thereunder the names and addresses of stockholders owning or holding one per cent 
or more of total amount of stock. If not owned by a corporation, the names and addresses of the 
individual owners must be given. If owned by a firm, company, or other unincorporated concern 
its name and address, as well as those of each individual member, must be given.) Owned by 
The Society of Exploration Geophysicists, Executive Committee: R. D. Wyckoff, Editor, Box 2038, 
Pittsburgh, Pa.; J. F. Gallie, Business Manager, Box 777, Austin Texas; W. T. Born, President, 
Box 2040, Tulsa, Oklahoma; H. B. Peacock, Vice-President, Houston, Texas; Andrew Gilmore, 
Secretary-Treasurer, Tulsa, Oklahoma; E. A. Eckhardt, Past-President, Box 2038, Pittsburgh, 
Pennsylvania. 

3. That the known bondholders, mortgagees, and other security holders owning or holding 
I per cent or more of total amount of bonds, mortgages, or other securities are: (If there are none, 
so state.) None; non-profit organization. 

4. That the two paragraphs next above, giving the names of the owners, stockholders, and 
security holders, if any, contain not only the list of stockholders and security holders as they appear 
upon the books of the company but also, in cases where the stockholder or security holder appears 
upon the books of the company as trustee or in any other fiduciary relation, the name of the person 
or corporation for whom such trustee is acting, is given; also that the said two paragraphs contain 
statements embracing affiant’s full knowledge and belief as to the circumstances and conditions 
under which stockholders and security holders who do not appear upon the books of the company 
as trustees, hold stock and securities in a capacity other than that of a bona fide owner; and this 
affiant has no reason to believe that any other person, association, or corporation has any interest 
direct or indirect in the said stock, bonds, or other securities than as so stated by him. 

5. That ithe average number of copies of each issue of this publication sold or distributed, through 
the mails or otherwise, to paid subscribers during the twelve months preceding the date shown 
above is not daily. (This information is required from daily publications only.) 

J. F. GALLIE 


}ss 


Sworn to and subscribed before me this 18th day of November 1940. 
[SEAL] ANNE W. GREEN 


(My commission expires June 1, 1941.) 





